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Abstract. We undertake a rigorous study of the stability of the Cauchy horizon 
in the naked self-similar Lemaitre-Tolman-Bondi spacetimes under even parity linear 
perturbations. We use a combination of energy methods and results about L p -spaces to 
determine the behaviour of the perturbations as they evolve through the spacetime. We 
first establish that an average of the perturbation generically diverges on the Cauchy 
horizon. We next introduce a rescaled version of the perturbation, and show that 
it is bounded and non-zero on the Cauchy horizon. This in turn shows that the 
perturbation itself diverges in a pointwise fashion on the Cauchy horizon. We give 
a physical interpretation of this result using the perturbed Weyl scalars. This result 
supports the hypothesis of cosmic censorship. 
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1. Introduction: Cosmic Censorship and Perturbation Theory 

The formation of naked singularities in various collapse models is a well known 
occurance; examples include the Reissner-Nordstrom (RN) and Kerr spacetimes pQ, 
certain models of critical collapse [2] and certain self-similar perfect fluid and dust 
solutions |3j. In response to the formation of such singularities, Roger Penrose proposed 
that the evolution of physically reasonable and generic initial data will not result in the 
formation of a naked singularity, a statement known as the cosmic censorship hypothesis 
(CCH) p]. Indeed, it is noted that many spacetimes which contain naked singularities 
also show high degrees of symmetry and it is tempting to conclude that the formation of 
the naked singularity is (in some cases) due to the unphysical symmetry of the spacetime. 

It follows that one manner in which the CCH can be studied is by disturbing 
the symmetry of a naked singularity spacetime by introducing perturbations. Should 
these perturbations remain finite as they evolve through the spacetime, then the 
singularity has displayed stability to such perturbations. In particular, the behaviour of 
perturbations on the Cauchy horizon (CH) of the RN spacetime is illustrative. Recall 
that the CH is the first null ray emitted by the singularity; it can equally be thought 
of as the boundary which separates observers who can see the singularity from those 
who cannot. In RN, the CH is a hypersurface on which perturbations diverge. More 
precisely, metric perturbations which arrive at the CH from the exterior have an infinite 
flux on the CH, as measured by observers crossing the CH. One might expect that 
perturbations evolving through other spacetimes containing naked singularities would 
show similar divergent behaviour on the CH. 

Should the perturbations turn out to behave in a finite manner on the CH, in certain 
cases one can still rule out the spacetime as a serious counter-example to the CCH. In 
some cases, the spacetime involves an unrealistic matter model, for example the Vaidya 
spacetime containing null dust or the perfect fluid spacetime, which neglects shear, 
viscosity and heat conduction. Another phenomenon which sometimes occurs is that 
the formation of the naked singularity depends on specific initial data; any perturbation 
of the initial data and the naked singularity fails to form [3]. 

We consider here the behaviour of even parity perturbations of the self-similar 
Lemaitre-Tolman-Bondi spacetime. We note that this spacetime cannot be taken as 
a serious counter-example to the CCH, as the first of the above-mentioned defects 
is present in this spacetime. The matter model used is a dust cloud which ignores 
pressure (and pressure gradients), and therefore does not provide a realistic description 
of gravitational collapse. Nonetheless, the simplicity of this spacetime makes it a very 
useful toy model for gravitational collapse resulting in naked singularity formation. 
For example, in the present paper we develop techniques that should in principle be 
applicable to any spherically symmetric self-similar spacetime. 

In a previous paper [5], the behaviour of odd parity perturbations of the self- 
similar LTB spacetime was considered. These perturbations were found to remain finite 
at the CH, where finiteness is measured with respect to certain energy norms of the 
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perturbation, and pointwise values thereof. This result holds for a general choice of 
initial data and initial data surface. In this paper, we use broadly similar techniques, 
with non-trivial extensions, to examine the behaviour of even parity perturbations of 
the same spacetime. 

The question of the behaviour of even parity perturbations of this spacetime was 
considered by Nolan and Waters j6]. In this work, both a harmonic decomposition 
and a Fourier mode decomposition of the gauge invariant perturbation variables and 
equations were used, and the individual Fourier modes X w ^ >m {z) were analysed. It was 
found that modes which are finite on the past null cone remain finite on the Cauchy 
horizon. However, the question of how to resum these modes on the Cauchy horizon was 
not fully addressed and so the problem of getting a complete analytic understanding 
the even parity linear perturbations remains unresolved. 

In the next section, we describe the structure of the self-similar LTB spacetime and 
examine conditions necessary for this spacetime to have a naked singularity. In Section 
[3J we discuss the perturbation formalism due to Gerlach and Sengupta. This formalism 
uses a multipole decomposition to write linear metric and matter perturbations in terms 
of functions depending on the similarity coordinate z and the comoving radial coordinate 
r. These perturbations do not share the background symmetries of self-similarity and 
spherical symmetry. In Section H] we derive a fundamental four-dimensional system of 
coupled linear PDEs which describe the evolution of the even parity perturbations. 

In practice, we work with a five dimensional system, which has the very useful 
property of symmetric hyperbolicity. These five variables obey an equation of motion 
arising from the linearised Einstein equations and must also satisify a constraint, the 
elimination of which produces the four variable system. Having determined a useful 
form for the perturbation variables and the Einstein equations, we aim to examine the 
behaviour of the perturbations as they reach the CH. In particular, we wish to know 
whether or not the perturbations remain finite there. 

Our strategy in addressing this problem is as follows. In Section HI we find that 
the perturbations obey a PDE whose coefficients, as a consequence of the self-similarity, 
depend only on the similarity coordinate. This means that we can take the state vector 
u which describes the perturbations and integrate with respect to the radius, to produce 
a kind of "average". Since the CH is a surface of constant similarity variable, and since 
the equations of motion have coefficients independent of the radius, we can reasonably 
expect that the behaviour of this averaged perturbation should reflect the behaviour of 
the perturbation itself. In other words, we use the behaviour of solutions to the averaged 
perturbation's equation of motion as a guide to the behaviour of the perturbation itself. 
This is the core strategy behind this work. 

In Section [5] we analyse the behaviour of solutions to the ODE which this averaged 
form satisfies, as they approach the CH. This ODE has a regular singular point, and 
methods for solving such systems are well understood. We can show that solutions 
to this system generically blow-up on the CH, with a characteristic power given by an 
eigenvalue of a particular matrix, which we denote — c. However, since this result applies 



Cosmic Censorship for Self-Similar Spherical Dust Collapse 



4 



only to an average of u, we cannot immediately conclude that u blows-up at the CH in 
a pointwise manner. 

In Section El we investigate the pointwise behaviour of u by introducing a new state 
vector x, related to u, which we expect to have a finite limit on the CH. We devote the 
next section to a series of results, which cumulatively show that x is finite and non-zero 
on the CH. This in turn establishes that u blows-up in a pointwise manner on the CH. 

Finally, in Section [7] we use the perturbed Weyl scalars to give a physical 
interpretation of our result. In Section [8] we make some concluding remarks and discuss 
further developments of this work. We use units in which G = c = 1. 



2. The Self-Similar LTB Spacetime 

2.1. The LTB Spacetime 

The Lemaitre-Tolman-Bondi spacetime is a spherically symmetric spacetime containing 
a pressure-free perfect fluid which undergoes an inhomogeneous collapse into a 
singularity. Under certain conditions this singularity can be naked. We will initially use 
comoving coordinates (t,r,9,4>), in which the radius r labels each successive shell in the 
collapsing dust. The line element for such a spacetime can be written as 

ds 2 = -dt 2 + e"(t, r)dr 2 + R 2 {t, r)dtt 2 , (1) 

where dQ 2 = d8 2 + sin 2 8d(p 2 and R(t,r) is the physical radius of the dust. The stress- 
energy tensor of the dust can be written as 

= p(t,r)W, 

where is the 4-velocity of the dust and p(t, r) is its rest mass density. In comoving 
coordinates, = Sq. 

The background Einstein equations for the metric and stress energy in comoving 
coordinates immediately provide the following results, 

„ /2 _ 1 dR 1 (dR\~ x dm 

6 ~ ^T+J(r) dr 7 P{ ,T) ~ AirR? \dr J dr ' () 

(£)"-^-/<o. (3) 

The function m(r) is known as the Misner-Sharp mass and is a suitable mass measure 
for spherically symmetric spacetimes. 

Recall that a shell focusing singularity is a singularity which occurs when the 
physical radius R(t, r) of the dust cloud decreases to zero, so that all the matter shells 
have been "focused" onto a single point. In this spacetime, a shell focusing singularity 
occurs on a surface of the form t = t s f(r), which includes the scaling origin (t, r) = (0, 0). 

We immediately specialise to the marginally bound case by setting f(r) = 0. See 
[5] for more details on the background spacetime. 
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2.2. Self- Similarity 

We follow here the conventions of |8J. A spacetime displays self-similarity if it admits a 
homothetic Killing vector field, that is, a vector field £ such that 

£|SV = 2 9w ( 4 ) 

In comoving coordinates, the homothetic Killing vector field is given by £ = t-jL + r4-. 
When self-similarity is imposed on the metric and stress-energy tensor, we find that 
functions appearing in the metric, the dust density and the Misner-Sharp mass have the 
following scaling behaviour 

v(t,r) = v(z), R(t,r)=rS(z), (5) 

^' r ) = _ ^2 - > m(r) = Xr, (6) 

where z — — t/r is the similarity variable and A is a constant (the case A = corresponds 
to flat spacetime). Having imposed self-similarity, we find that combining (j3J), flSJ) and 
(JE]), produces 

S(z) = {az + 1) 2/3 , (7) 

where a = 3y|. The form of the line element (pQ) is invariant under the coordinate 
transformation r — > r — f(r). We choose r so that R\t=o — r and |jjr|t=o = 1- With this 
expression for S(z) we can solve for ^ explicitly In (j2J) we convert ^ to a derivative 
in (z,r) and find that 

e^ = ^ = ( 1 -az + l)(l + az)- 1 / 3 . (8) 

We note that by combining §5§ with (j2j), we can find an expression for q(z), 

a 2 

q ^ Z) = 6tt(3 + Aaz + o? z 2 )' (9) 
We state the metric in (z, r) coordinates for future use, 

ds 2 = -r 2 dz 2 + e u (z)(l - z' ! e~ l '(z))dr 2 - Irzdrdz + R 2 dVi 2 . (10) 

In Section [7] we will need the radial null directions of the self-similar LTB spacetime. In 
terms of (z, r) coordinates, retarded null coordinates u and v take the form 

" = rexp (-f /So)' " = rexp ("f 7§))' 

where f± := ±e u ^ 2 + z. In these coordinates, the metric takes the form 

t 2 

ds 2 = (1 -e v z- 2 )dudv + R 2 (t,r)dn 2 . 

uv 

In order to calculate the perturbed Weyl scalars, we will need the radially in- and 
outgoing null vectors, Z M and n M . These vectors obey the normalization g liu l^"n u — — 1. 
A suitable choice is therefore 

B{u,v)ou ov 
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where B(u,v) = ^ ^1 — s -^rj- In what follows, we shall always use a dot to indicate 
differentiation with respect to the similarity variable z, ■ = J^. 

2. 3. Nakedness of the Singular Origin 

We now consider the conditions required for the singularity at the scaling origin 
(t, r) = (0, 0) to be naked. As a necessary and sufficient condition for nakedness, the 
spacetime must admit causal curves which have their past endpoint on the singularity. 
It can be shown [S] that it is actually sufficient to consider only null geodesies with 
their past endpoints on the singularity, and without loss of generality, we restrict our 
attention to the case of radial null geodesies (RNGs). 

The equation which governs RNGs can be read off the metric (pQ), 

% = ±e^ 2 . (13) 
dr 

Since we wish to consider outgoing RNGs we select the + sign. We can convert the 
above equation into an ODE in the similarity variable, 

z + rz' = -e v ' 2 . (14) 

We look for constant solutions to this equation, which correspond to null geodesies that 
originate from the singularity. In other words, the existence of constant solutions to 
(fT4l) is equivalent to the nakedness of the singularity. For constant solutions, we set the 
derivative of z to zero, so that the Cauchy horizon is the first root of 

z + e^ 2 = z + S{z) -zS{z) =0, (15) 

where we reserve the dot notation to indicate a derivative with respect to z, • := d/dz. 
By using (J7J), we can write (jl~5l) as an algebraic equation in z, 

az" + (\ + ^jz 3 + (^j z 2 + az + 1 = 0. (16) 

We wish to discover when this equation will have real solutions. This can easily be 
found using the polynomial discriminant for a quartic equation, which is negative when 
there are two real roots. In this case we have 

D = 4=(-729 + 2808a 3 -4a 6 ), (17) 

— I 

which is negative in the region a < a* where a* is 

3 , . 

0.638... 18 



(2(26 + 15VS)) 1 / 3 

This translates to the bound A < 0.09. From (|6]), we can see that this result implies 
that singularities which are "not too massive" can be naked. See Figure [1] for a Penrose 
diagram of this spacetime. 



Remark 2.1 In fact, one can find D < in two ranges, namely a < a* « 0.64 and 
a > a** ~ 8.89. We reject the latter range as begin unphysical. Consider ([7]), which 
indicates that the shell-focusing singularity occurs at z — — 1/a. If we chose the range 
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Figure 1. Structure of the self-similar LTB spacetime. We present here the conformal 
diagram for the self-similar LTB spacetime. The gray shaded region represents the 
interior of the collapsing dust cloud. We label the past null cone of the naked singularity 
by J\f, future and past null infinity by J + and J~ and the Cauchy horizon by H. 

a > a** we would find that the corresponding outgoing RNG occurs after the shell 
focusing singularity and so is not part of the spacetime. 

Remark 2.2 We note that this analysis has assumed that the entire spacetime is filled 
with a dust fluid. A more realistic model would involve introducing a cutoff at some 
radius r = r*, after which the spacetime would be empty. We would then match the 
interior matter-filled region to an exterior Schwarszchild spacetime. However, it can be 
shown that this cutoff spacetime will be globally naked so long as the cutoff radius is 
chosen to be sufficiently small [7]. We will therefore neglect to introduce such a cutoff. 

3. The Gerlach-Sengupta Formalism 

We shall use the Gerlach-Sengupta method [9] to perturb this spacetime (we follow 
the presentation of [TO])- This method exploits the spherical symmetry of the 
spacetime by performing a decomposition of the spacetime into two submanifolds (with 
corresponding metrics). Perturbations of the spacetime are then expanded in a multipole 
decomposition and gauge invariant combinations of the perturbations are constructed. 
We begin by writing the metric of the entire spacetime (J^A 4 ,g^ u ) as 

ds 2 = g AB (x c )dx A dx B + R 2 (x c )>y ab dx a dx b , (19) 
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where qab is a Lorentzian metric on the 2- dimensional manifold M 2 and 7 a b is the metric 
for the 2-sphere S 2 (and the full manifold is M 4 = M 2 x S 2 ). The indices A,B,C... 
indicate coordinates on A4 2 and take the values A, B . . . = 1, 2 while the indices a,b,c. . . 
indicate coordinates on S 2 and take the values a, b . . . = 3, 4. The covariant derivatives 
on Ai A , Ai 2 and S 2 are denoted by a semi-colon, a vertical bar and a colon respectively. 
The stress-energy can be split in a similar fashion, 

t^dx^dx" = t AB dx A dx B + Q(x c )R 2 'y ab dx a dx b , 

where Q(x c ) = \t a a is the trace across the stress-energy on S 2 , which vanishes in the 
LTB case. Now if we define 
R\A 



V = -^ + 2v A A + 3v A v A , 

then the Einstein equations for the background metric and stress-energy read 
Gab = -2(v a \b + v A v B ) + VqQab = 87it AB , 

l -G\ = -H + v A v A + v A A = 8ttQ(x c ), 

where G a a = 7 ab G a 6 and TZ is the Gaussian curvature of Ai 2 , TZ = \R ( ~a >A where R^ B 
indicates the Ricci tensor on M 2 . 

We now wish to perturb the metric (fT9|) . so that g^ u (x a ) —> g^ u (x a ) + 5g^ u {x a ). To 
do this, we decompose Sg^ u (x a ) and write explicitly the angular dependence using the 
spherical harmonics. We write the spherical harmonics as Y™ = Y. {Y} forms a basis 
for scalar harmonics, while {Y a := Y a ,S a := e b a Yb} form a basis for vector harmonics. 
Finally, {Y^ a b, Z ab := Y a ,b + ^^ Y^ab, S a -b + Sb- a } form a basis for tensor harmonics. 

We can classify these harmonics according to their behaviour under spatial inversion 
x — x: A harmonic with index I is even if it transforms as (—1)' and odd if it 
transforms as (— According to this classification, Y, Y a and Z a b are even, while 
S a and S( a: b) are odd. 

We now expand the metric perturbation in terms of the spherical harmonics. Each 
perturbation is labelled by (/, m) and the full perturbation is given by a sum over all I 
and m. However, since each individual perturbation decouples in what follows, we can 
neglect the labels and summation symbols. The metric perturbation is given by 

SgAB = h A BY, 5g A b = h A Y :b + h%S b , 

8g ab = R 2 KY lab + R 2 GZ ab + h(S a:b + S ha ), 

where Hab is a symmetric rank 2 tensor, h A and h® are vectors and K, G and h are 
scalars, all on M 2 . We similarly perturb the stress-energy t^ u — > t^ u + 5t fll/ and expand 
the perturbation in terms of the spherical harmonics, 

5t AB = At AB y, 5t Ab = At*iY : b + At%S b , (20) 
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5t ab = r 2 At 3 lab Y + r 2 At 2 Z ab + 2AtS {a:b) , (21) 

where At AB is a symmetric rank 2 tensor, At^ and At° are vectors and At 3 , At 2 and 
At are scalars, all on M 2 . 

We wish to work with gauge invariant variables, which can be constructed as 
follows. Suppose the vector field £ generates an infinitesimal coordinate transformation 
x — > x' = £ + We wish our variables to be invariant under such a transformation. We 
can decompose £ into even and odd harmonics and write the one-form fields 

f = U{x c )Ydx A + e(x C )Y-.adx\ f = £°S a dx a . 

We then construct the transformed perturbations after this coordinate transformation 
and look for combinations of perturbations which are independent of £, and therefore 
gauge invariant. We will list here only the even parity gauge independent perturbations 
and we set the odd parity perturbations to zero. In the even parity case, the metric 
perturbation is described by a gauge invariant 2-tensor k AB and a gauge invariant scalar 

K 

kAB = h AB - (p A \ B + pb\a), k = K - 2v A p A , (22) 

where p A = h\ — (1/2)R 2 G\ A . The even parity gauge invariant matter perturbation is 
given by 

Tab = At AB - t AB lC Pc ~ 2(t CA p c w + t C BP C \ A ), (23) 



T A = At A - t A c pc - R 2 (^j G\ A , 



(24) 



At 3 - t- ( , T 2 = At 2 — ( ^ ) G. (25) 



r 2 v 2 ; |c ' ~ v 2 

In the Regge- Wheeler gauge, h B — h — G — 0, which implies that pa = and the 
gauge invariant matter perturbations coincide with the bare matter perturbations, which 
simplifies matters considerably. 

We now list the linearised Einstein equations, which can be written in terms of the 
gauge invariant quantities listed above. We find 

C „ /<->;„ \D , n \C 



\o v 

' AB 



(k C A\B + k C B\A - k A B\c)v ~ g AB (2k CD l ~ k D 

~{k\AVB + k\ B v A + k\ AB ) + (vq + ^2ffi^ ) k/ 

-9ab (V + lk DF v D \ F + 3k DF v D v F ^j 

+9AB ( {l ~ 1 ^2 +2) k - k \F - 2 V« F ) = to^AE, 

- k AB lAB + (*a)\b - lk AB \ A v B + k\ B v B + R AB (k AB - kg AB ) 

--^r- kA A + k \A+ U lA v A = 16vrT 3 , (27) 

kJ B -k B B \ A + k B B v B - k\ A = IQttT a , (28) 
k A A = 16ttT 2 . (29) 



(26) 
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Finally, we present the linearised stress-energy conservation equations, for completeness. 
(t AB \ D + 2t AB v D )k BD + Q{k\ A - 2kv A ) - t AB k lB + h BC k BC]A - h AB kp B 

+t AB k BF lF , -j- 2 (R 2 T AB )\ B + l -^^T A + 2v A T 3 = (30) 

^T B r + T 3 — (/ " 1} 2 (/ + 2) g = \k AB t AB + Q(* - ^/). (31) 
We now consider the form of these equations for the background LTB spacetime. 

3.1. The Matter Perturbation 

To proceed further, we must find a relation between the gauge invariant matter 
perturbation terms (1231 - [25]) and the dust density and velocity discussed in section 
2.1. This amounts to specifying the matter content of the perturbed spacetime. To 
do this, we write the stress-energy of the full spacetime as a sum of the background 
stress-energy and the perturbation stress-energy (where a bar indicates a background 
quantity), 

We assume that the perturbed spacetime also contains dust and write the density as 
p = p + 5p and the fluid velocity as w M = + 5u^. We can now find an expression for 
the perturbation stress-energy (keeping only first order terms), 

5T^ V = ~p(u ll 5u v + Uyb'u^) + b~pu^u v . (32) 

The perturbation of the dust velocity can now be expanded in terms of spherical 
harmonics as 5u^ = (5u A ,5u a ) = (5u A Y,SuEY :a ). By imposing conservation of stress- 
energy and requiring that the perturbed velocity = u fl + Su^ obeys u^u^ — — 1, one 
can show that the perturbed velocity can be written in the form 

<^=(<^r(t,r)y,7(t,r)K a ), (33) 

where the variable V acts as a velocity potential and obeys an equation of motion arising 
from perturbed stress-energy conservation (specifically, from the acceleration equations 
arising from stress-energy conservation), 

dT 1 

S = - 5 «(»,r), (34) 

where we have labelled the first component of k AB (see ff22l ) as k 00 = a(z,p). In 
addition, by using (I3T1) one can show that 

r y(z,p)=T(z,p)+g(p), (35) 

where p = lnr and g(p) is an initial data function for the velocity perturbation. If we 
compare this form for the perturbed stress-energy to the Gerlach-Sengupta form, we can 
find the gauge invariant matter perturbations for the LTB spacetime, which we write in 
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terms of (z,p) coordinates. The gauge invariant tensor Tab is given in (z,p) coordinates 
by 

T Q0 = 2pe p — + e 2p 5p, 

T 01 = T 10 = e p P (^ + z^) + e 2p z6 P , 

or 

T 11 = 2pze p — + e 2p z 2 Sp. 
op 

The vector T4 is given by Ta = (e^ p r y(z, p) , e p p r y(z, p)z) and the gauge invariant scalars 
both vanish, T 2 = 0, T 3 = 0. 



4. Reduction of the Perturbation Equations and the Main Existence 
Theorem 

In what follows, we omit the exact form of various matrices and vectors if they appear 
in versions of the system of perturbation equations which we do not use; relevant terms 
are included in Appendix A as indicated. As we have imposed self-similarity on this 
spacetime, a natural set of coordinates to present the linearised Einstein equations in is 
x^ 1 = (z,p, 8, </>), where z = — t/r is the similarity coordinate and p = ln(r). The metric 
in these coordinates takes the form 

ds 2 = e 2p {-dz 2 + {{S - zS) 2 - z 2 )dp 2 - Izdzdp + S 2 dQ 2 ). 

Our initial full set of perturbation equations consists of both components of ( 1271) . the 
t—p and p — p components of ( 126]) and the equation of motion (l34"j) for T(z,p). We now 
discuss the series of simplifications which allows us to make this choice, before stating 
the initial six dimensional system of mixed evolution and constraint equations 0- 

(1) We note that since T 2 = in this spacetime, by ( 129]) the metric perturbation tensor 
kAB is trace-free. We use this to eliminate one component of Uab- 

(2) Since T 2 vanishes, ( J29l) implies that the metric perturbation kAB is trace-free. 
Additionally, given that the scalars Q and and T 3 both vanish in the LTB case, one 
can show that (J27J is identically satisfied, assuming that the background Einstein 
equations, ( 1281) and ( 13T1) all hold. We will therefore use (}3T|) in preference to ( |27j) . 

(3) We note that the t — t component of ( 1261) gives us a relation for the perturbation of 
the dust velocity 5p in terms of the velocity and metric perturbation. We use this 
equation to eliminate 5p from the system. 

(4) The variable k is the only variable which appears at second order in derivatives 
in the resulting system. We therefore introduce a first order reduction by letting 

u(z,p) = k(z,p), v(z,p) = dk/dp and w(z,p) = k, where • = d/dz. 

f To download a Mathematica notebook containing these calculations, go to 
www.student.dcu.ie/~duffye27. 
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(5) The resulting system of equations consists of both components of (KT]) . the t — p 
and p — p components of (126]) and the equation of motion flMl) for T(z,p). To this 
we append the auxiliary equations k = w and v = dw/dp, which makes a total of 
seven equations. 

(6) This results in a first order system of seven equations for six variables, which we 
combine into a vector 

X = (a(z,p),/3(z,p),u(z,p),v(z,p),w(z,p),T(z,p)) T G R 6 . 

Here a(z,p) and f3(z,p) are the independent components of the tensor k^B i n {z,p) 
coordinates 



kAB 



a(z,p) /3{z,p) 



and since T 2 = 0, ([29]) implies that 5 = 2zf3(z,p) + (S(S -2zS) + z 2 (-l + S 2 ))a(z,p). 
Finally T(z,p) is the velocity potential given in ( 133]) . 

(7) We can write the system of equations in a more compact form as 

f) X f)X 

M{z)°— + N(z)^- + 0(z)X = S 7 (z,p), (36) 

for 7x6 matrices M(z), N(z) and O(z) and a 7 dimensional source vector 5*7. 
The dimensions of this system suggest that it may be possible to rewrite it as a 
6 dimensional system of evolution equations with a constraint. To identify the 
constraint, we look for linear combinations of the rows of M(z) which add to give 
zero. This corresponds to a linear combination of the equations in (136]) which has 
no time derivatives, that is, a constraint. We call this the Einstein constraint. 

(8) Having identified this constraint, we construct a new system consisting of six of the 
equations of the original system, with the constraint added to them. This is fully 
equivalent to the original seven equation system. We can write this system in a 
similar manner to (1361). 



f)X f)X 

P{z)— + Q(z)— + R(z)X = E 6 (z,p), (37) 

for 6x6 matrices P{z), Q(z) and R{z) and a source vector Eg. We note that in 
addition to the Einstein constraint, we have the trivial constraint dk/dp = v(z,p). 

(9) Having identified the two constraints, we would like to use them to eliminate two 
variables from the system and thus reduce the number of variables from six to four. 
To do this, we need to diagonalise the system. We multiply through (137]) by P" 1 and 
find the Jordan canonical form T of the matrix coefficient of dX/dp, T := P~ l Q. 
We also identify the similarity matrix S such that T = S^^^TS. 

(10) To diagonalise the system, we let Y = SX, where det(S') ^ 0. Y obeys the equation 
dY - dY 

— + T{z)— + R{z)Y = a 6 (z,p), (38) 
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where T = S^TS, R = S + S^P^RS and cr 6 = S' _1 P _1 S 6 . The matrix T is in 
Jordan form but is not diagonal. 

'11) In terms of the components of X, the components of Y G M 6 are 

Yi{z,p) = —a(z,p) - 8nzq(z)T(z,p) + v(z,p) 

S(z) 

-zw(z,p) + I^t/^p)' 

Y 2 (z,p) = 8Tcq(z)T(z,p), 
Y 3 (z,p) = u(z,p), 

-f3(z,p) + S(8irq(z)T(z,p) + w(z,p))(S - zS) 



Y 4 (z,p) 



S(S — zS) 

[z + S(S — zS))a(z,p) 



Y 5 (z,p) 



Y 6 (z,p) 



S(S — zS) 
(1 + S)(P(z,p) + (-z-S + zS)a(z,p)) 

2S{S - zS) 
(-1 + S)(-P(z,p) + (z-S + zS)a(z,p)) 
2S{S-zS) 

where S(z) is the radial function and q(z) is the density given by (jHJ). 

(12) The Einstein constraint can be given in terms of Y as 

Cl (z)Yx(z,p) + c 2 (z)Y 2 (z,p) + c 3 (z)Y 3 (z,p) + c 4 (z)Y 4 (z,p) (39) 

dY 

+c 5 (z)Y 5 (z,p) + c G (z)Y 6 (z,p) + c 7 (z) — {z,p) + c 8 (z)g(p) = 0. 

where the coefficients are not important. The g(p) term comes from the source term 
Tjq in (1361) . which could be written as = b(z)g(p) where the form of b(z) is not 
needed. 

(13) The trivial constraint dk/dp = v(z,p) can be stated in terms of Y as 

Yfap) + zY 4 (z,p) + (z-S + zS^j Y 5 (z,p) (40) 

/ \ BY 

+ (z + S-zS) Y 6 (z,p) - -g^(z,p) = 0. 

(14) Lemma 4.1 Suppose that f3ty) is satisfied on an initial surface z = zq. Then 
assuming that ( tffflj] holds, ^39$ will be satisfied on all surfaces z G (z c ,Zq\. 

Proof This is a straightforward but lengthy calculation which was carried out using 
computer algebra. □ 



This lemma indicates that the constraint ( |39l) is propagated by the system. 

(15) The trivial constraint (|4"0"1) is also propagated in the sense of Lemma 14. 11 The 
existence of these two constraints suggests that the true number of free variables 
in this system is four. We therefore aim to reduce this system to a free evolution 
system of four variables by using the two constraints to eliminate two variables. 
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This is carried out in such a way as to keep the system always first order in all 
variables. 

(16) We carry out this reduction in two steps, by first passing to a five dimensional system 
(which is still a mixed evolution-constraint system) by solving the trivial constraint 
for one variable, and then passing to a four dimensional system by solving the non- 
trivial constraint for another variable. The advantage to carrying out the reduction 
in this manner is that the five dimensional system is symmetric hyperbolic, which 
will be useful in what follows. 

(17) The five dimensional system is obtained by solving the trivial constraint ( 140]) for the 
variable Yi(z,p). We eliminate this variable and reduce the system to five variables, 
with a state vector w G M 5 . Again we would like to put the system in a diagonalised 
form. We do this as before by calculating the Jordan canonical form of the matrix 
coefficient H(z) of dw/dp and letting u = Sw e R 5 , where S is the similarity matrix 
arising from the transformation of H{z) into Jordan canonical form. In terms of Y, 
the new variables are given by 

ui(z,p) = f 1 (z)Y 3 (z,p) } 

u 3 (z,p) = Y A (z,p) + f 2 (z)Y 3 (z,p), 

u s (z,p) = Y 6 (z,p) + f 4 (z)Y 3 (z,p), 



u 2 {z,p) = Y 2 (z,p), 
u 4 (z,p) = Y 5 (z,p) + f 3 (z)Y 3 (z,p), 



where 



h{z) 



hi?) 



1 - S(z) 
2S(z) 

1 + S(z) 
2S(z) 



-zS{zf + S(z)(S(z) + zS{z)) 
S{z){S{z) - zS(z)) ' 
-1 + S(z) 



Hz) 



2S(z) 



The five dimensional system obeys an equation of motion of the form 

E 5 (z,p). 



_ + A{z) _ + C{z )u 



(41) 



We shall use this form of the system in the next two sections chiefly because it has 
the useful property that it is symmetric hyperbolic. This means that the matrix A, 



A 



/ 























/ 







-(z-S + zS)- 1 

-(z + 3- zS)~ z 

is symmetric. The matrix C and the source term S5 are given in Appendix A. 
(19) In terms of these new variables, the Einstein constraint becomes 

g 1 {z)u 1 {z,p) + 92(z)u 2 (z,p) + g 3 (z)u 3 (z,p) + g 4 (z)u 4 (z,p) (42) 

+g 5 (z)u 5 (z,p) + g 6 (z)-^-(z,p) + g 7 (z)g(p) = 0, 

where the coefficients gi{z), i = 1, . . . , 7 are listed in Appendix A, and g(p) is an 
initial data function. We note that g${z) vanishes on the Cauchy horizon. 
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(20) In order to eliminate one more variable, we solve f !42p for w 2 (z,p). As before, we 
then put the new system in Jordan canonical form by writing it in terms of the 
vector k E M 4 , where in terms of u(z,p), the new variables are 

u 3 (z,p) 



where 

/(*) = 



h(z,p) = ux(z,p), k 2 (z,p) , . 

/(*) 

k 3 (z,p) = u 4 (z,p), h(z,p) = u 5 (z,p), 



12(3 + 2az)(2a - 3h(z)) 2 (2a + 3h(z)) 



(3 + az)(16a 4 + 108a 2 z - 81h(z) - A8a 3 h(z) - 27a(-4 + 3zh{z))) 
and h(z) = (1 + az) 1/3 . 
(21) k obeys the differential equation 

Cfn (in — * — * 

— + E(z)— + B(z)k = i: 4 (z } p) J (43) 

— * 

where E{z) is given in Appendix A and we omit B(z) and £4. This system is a 
free evolution system in the sense that there are no further constraints which must 
be obeyed by these variables. The system cannot be reduced to any simpler form 
than this. However, the matrix E(z) is not symmetrizable, which implies that this 
system is not symmetric hyperbolic. This is why we choose to work with the five 
dimensional system ( T4Tj) and the Einstein constraint (f4"2]) . 

We will slightly rewrite the five dimensional system (jH]) as 

(111 f)i7 -, 

t— + A{t)— + C{t)u = nt,v) (44) 

where now t = z — z c , so that t = is the Cauchy horizon. In terms of the coefficient 
matrices and source in (llTi) . A(t) = tA(z), C(t) = tC(z) and £(£,£>) = tYl^{z,p). We 
briefly list here the most important properties of this system. 

• u(t,p) is a five dimensional vector, whose components are linear combinations of 
the components of the gauge invariant metric and matter perturbations. 

• A(t) and C(t) are five-by-five matrices. We note that A(t) = tA(z) and the 
matrix A(z) contains a factor of h~ 1 (z) in the (5,5) component. Here h(z) : = 
z+S — zS = z—f(z), where S(z) is the radial function (see (JTj)) and f(z) = —S+zS. 
f(z c ) = z c so that h(z) vanishes on the Cauchy horizon. If we Taylor expand 
h(z) = z — z c — f(z c )(z — z c ) + 0((z — z c ) 2 ) = t — f(z c )t + 0(t 2 ), then we can see 
that th~~ 1 (z) is analytic at the Cauchy horizon where t = 0. This in turn implies 
that A(t) = tA(z) is analytic at the Cauchy horizon. Similar remarks apply to the 
matrix C(t) = tC(z), since the fifth row of C(z) contains h~ 1 {z) factors. Similarly, 
the fifth component of £5 contains a factor of h^ 1 {z). So overall, A(t), C(t) and 
£(i,p) are analytic for t > 0. 

• A(t) is diagonal, whereas C(t) is not. The first four rows of A(t), C(t) and £(t,p) 
are 0(t) as t — > 0, while the last row of each is 0(1) as t — > 0. 
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• The source S is separable and we can write it as H(z,p) = h(t)g(p), where h(t) is 
an analytic vector valued function of t and g(p) is an initial data term. 

• g(p) represents the perturbation of the dust velocity and is a free initial data 
function. The results below follow for g 6 Cq°(M., R) which we assume henceforth. 

In what follows, t\ is the initial data surface and < t < t%, where t = is the Cauchy 
horizon, so the Cauchy horizon is approached in the direction of decreasing t. 

Theorem 4.2 The IVP consisting of the system fij4\ ) along with the initial data 

u = f(p) 

H 

where / 6 C£°(R, 1R 5 ), possesses a unique solution u(t,p), u e C°°(R x (0,ti],R 5 ). For 
all t e (0, t\], u(t, •) : R — > R 5 has compact support. 

Proof This is a standard result from the theory of symmetric hyperbolic systems, see 
Chapter 12 of [II]. □ 

We note that since the constraint is propagated by the five dimensional system (see 
Lemma ETTi) . a choice of smooth and compactly supported initial data for the components 
Ui(z,p), u 3 (z,p), u^Zip) and u 5 (z,p) is sufficient to ensure that u 2 (z,p) as given by the 
constraint ()42l) is also smooth and compactly supported. Therefore, this theorem also 
provides sufficient conditions for the existence of unique solutions to the four dimensional 
free evolution system. 

Remark 4.1 In Section [6] we will require solutions u(t,p) in L 1 (R, R 5 ) for a choice of 
initial data / G L 1 (R, R 5 ). It follows immediately from Theorem 14.21 by the density 
of in L 1 that for < t < t u u{-,p) e L X (R, R 5 ). To show that we can extend 
our choice of initial data to L 1 , we require a bound on u, which is established in the 
following lemma. For this lemma, we will need Gronwall's inequality [12], which states 
that for continuous functions <f)(t), ip(t) and xW? if 

<V(*)+ / X(s)<p(s)ds, 

J a 

then 

(pit) < 4>{t) + J x(s)ip(s)exp (J^ x(v>)du\ ds. 

Lemma 4.3 The I}-norm ofu(t,p) obeys the bound 

IKOIIi^cx^lKtOHx + ca^lHl!, (45) 
for < t < ti, where C\(t) and 02(f) are continuous on (0, ti]. 



Proof u obeys the equation 

du . , . du 
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where £ = h(t)g(p) and g(p) in an initial data term. Each row of (Ho]) can be written 

as 

t-^- + a^t)-^ 1 + Ci(t)ui = T,i(t,p) - 22 °i3 ( ( t ) u A t ^P) = Si(t,p), (47) 

for i = 1, . . . , 5, where eii(i) and Cj(t) are the diagonal components of the matrices A(t) 
and C(t) respectively, and because C(t) is not diagonal, the off-diagonal components 
Cij(t) are put into the source term Si(t,p). We can solve (1471) using the method of 
characteristics. The characteristics are given by 

^ = ^ =► Pi (t) = Vi + n i (t), (48) 
at t 

where 7Tj(t) = — f* 1 ^^-dr and = Pi(ti). On characteristics, ( I4"7j) becomes 

t-^(t, Pi (t)) + c i (f)u ! (t,p I (t)) = Si(t, Pi (t)). (49) 

The integrating factor for fT49|) is e^(i) where £j(t) = — f* 1 ^p-dr, and the solution to 
(1491 is 

Wi^Pi) = e-^(t) M S 0) ( Pi - 7Ti(*)) - e 6 (*) / 1 l^l^^pi + TTj (t) - 7ri(t))dr. (50) 

We take the L 1 -norm by taking an absolute value and integrating with respect to p; this 
produces 

+e^(t) jf 1 ^M|^(r,p, + n{r) - tt^))^ dp. 

Recall Theorem O which tells us that at each * G (0,*i], m,- G Cg°(M,M). This allows 
us to apply Fubini's theorem; that is, to interchange the order of the integrals above. 
We note that the structure of the characteristics (|4"8"1) indicates that evaluation of the 
L 1 -norm of a function f(t,p) at fixed time t yields the same result as the evaluation of 
the L 1 -norm of / evaluated on characteristics, that is 



\f(r,p\dp= / \f(r,Pi + 7T;(t) - 7T»(t))| dp. 
So, applying Fubini's theorem and using the form of Si produces 

II^WIIi = e^(t)IK (0) ||i+ (51) 
<*(f) (W)IIMIi+ E Mr)\ [ \u J \dp)dr. 

V 3=1,3& JR / 



T 

Then if we note that \\g\\i does not depend on t, we can write 

rh 

IK^Hx = sup <d 1 {t)\\^\\ 1 + d 2 {t)\\g\\ l + d 3 {t) / d^Wuir^dr, (52) 



Cosmic Censorship for Self-Similar Spherical Dust Collapse 



18 




Figure 2. Characteristics: We show here typical characteristic curves pi — Pi(t), along 
which the solution ([50)) is evaluated. 



where the di(t) functions are the suprema of the various t-dependent functions which 
appear in (15 ip (and their precise value is not important). Now applying Gronwall's 
inequality (with ip(t) = di(t)\\u^\\i + c^WIMIi) produces 

\\^\\i<d 1 (t)\\u^\\ 1 + d 2 (t)\\g\\i 

+ J t 1 d 5 (r) (d 1 (r)\\u^\\ 1 + d 2 (r)\\g\\ 1 ) exp U * d 5 {r')dA dr, 

where the exact value of d^{t) is unimportant. Again, since ||g||i and ||{r^||i do not 
depend on t, we can summarize this as 

IK*)||i <ci(t)|Kti)||i + c 2 (t)||^|| b 

for some continuous functions Ci(t) and C2(t) defined on (0, ti], whose exact value is not 
important. 

□ 



Remark 4.2 The characteristics Pi(t) provide a C 1 foliation of the region Q = {(t,pi) : 
t E (0, ti],Pi E M}. A typical such foliation is shown in Figure |2j For every point 
q = {t2iPi,2) E Q-i there is a unique characteristic C\ such that q E C\ which we label 
by r]i,2 = Pi(h) • Define the set Q t '■= {(t',p'i) E £1 : f = i}. Then the characteristics 

Ci 

provide a natural diffeomorphism of fit, 
Pi,2 E fif 2 -> Pi,3 E fl t3 , 

with 

Pi,2 = Pi,3 ~ / ai(r)dr. 
Jt 2 

The fact that u obeys a bound of the form (l4"5j) allows us to extend our initial data to 
f = u(t 1 ,p)EL 1 (M 
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Theorem 4.4 The IVP consisting of the system along with the initial data 
u = f(p) 



where / G L^IRjIR 5 ), possesses a unique solution u(t,p), u G C°°((0, tx], L l (R, R 5 )). 

Proof The proof of this result relies on the bound ( 1451) and follows by a standard 
argument exploiting the density of C^°(R, 1R 5 ) in M 5 ). See Theorems 5 and 7 of 

[13] for examples of such techniques and Chapter 12 of [XT] for background details. □ 

Although these theorems provide for the existence of smooth or L 1 solutions prior to 
the Cauchy horizon, it gives us no information about their behaviour as they reach 
the horizon itself. We must therefore consider this behaviour separately. We note that 
this problem is rendered nontrivial by the fact that the Cauchy horizon is a singular 
hypersurface of (1141) . 

5. Behaviour of the L 9 -Norm 

In this section, we will use Theorem 14.41 to provide for the existence and uniqueness of L 1 
solutions to ( jUl) with a choice of L l initial data. This theorem applies for t G (0, tx] only 
and we must consider the behaviour of u on the Cauchy horizon separately. Our strategy 
in tackling this problem is as follows (see [H]). We expect that any divergence which 
might arise in the perturbation would be in some sense (to be defined) independent 
of the radial coordinate, since the Cauchy horizon is a hypersurface of constant t, and 
since the coefficients of ( 1H1) are independent of p. Motivated by this observation, we 
introduce the integral of the perturbation vector with respect to the radial coordinate, 
which acts as a kind of "average" of the perturbation. This variable obeys a relatively 
simple system of ODEs, the solutions to which can be determined. 
Let u be the solution of (jH|) with u(tx,p) = u^ '- Then define 



which is a kind of "average" of u(t,p) (note that the existence of u is guaranteed since 
M < ||w||i < oo since u G L 1 ). If we integrate with respect to p through the system 
( HID , we find that u obeys the ODE 



(see Chapter 2 of [15] and Chapter 4 of [12] for the theory of such points). We now 
state a theorem which gives the fundamental matrix for this system. Recall that the 
fundamental matrix for an ODE system is a matrix whose rows are linearly independent 
solutions to the ODE in question. 




(53) 
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Theorem 5.1 The fundamental matrix corresponding to (5J\) is 

H{t) = J(t) + K[t) (55) 

= p{t)r Co +p(t)r Co / p-^rjT^-^dr, 



where P(t) = I+tP\ + . . . is a matrix series whose coefficients can be found by a recursion 
relation from the Taylor expansion ofC(t). Cq is the Jordan canonical form of the zero 
order term in the Taylor expansion of C(t). Co takes the form Cq = diag(0,0,0,0,c) 
where c is a constant given by 



c = lim t 



3 + {S{z c ) - z c )S{z c ) + S{z c ){-3 + z c S{z c )) 



^° V {l-S{z c )){z c + S{z c )-z c S{z c )) 
c depends only on a, and satisfies c G (3, +oo) ; with 
lim c = 3, lim c = +oo, 

where a* is the maximum value of a for which a naked singularity forms. 



Proof The proof is a standard result for systems of the form flMl) . See [15] or [12] for 
details. We have relegated the proof of the results about the behaviour of c to Appendix 
B.Q □ 

Our next task is to analyse in more detail the behaviour of this fundamental matrix. 



5.1. Behaviour of the Fundamental Matrix 

We now expand about the Cauchy horizon so that P(t) = I + 0(t) and write the 
homogeneous part and the particular part of H(t) separately. Given the form of the 
matrix Co listed in Theorem 15.11 the homogeneous part takes the form 

J(t) = diag(l + o(t), 1 + o(t), 1 + o(t), i + o(t), r c + o(r c+1 )), 

and the particular part can be written 

K(t) = diagOi, k 2 , «3, «4, «b), 

where 

rti 

Ki = J T- X %(T){l + 0{T))dT, 

for i — 1, 2, 3, 4 and 

k 5 = t~ c J 1 r^ 1 E 5 (r)(l + 0{t)) dr. 

f We note that there exists a set A of values of a such that c(a) is a natural number, that is 
A = {a e (0,o*) : c(a) £ Nfl (3,oo)}. When c e Nfl (3,oo), the fundamental matrix ([55]) will 
contain extra log terms. However, since this set has zero measure in the set a £ (0,a*) we will not 
consider it further. See [15] and [12] for further details. 
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We will now examine the particular part. We recall that £, is separable, so that 
£ = h(t)G, where h(t) = (0, 0, tk^(t), tk&(t), ks(t)). The fej(i) functions are all analytic 
at t = and G := j R g(p) dp G R. Now the k terms become 

Ki = /t 2 = 0, 

Kj = Gkj(t*){ti - t) + 0((ti - t) 2 ), 
for j = 3, 4 and 

K5 = r c k 5 (t*)G f^-^J + o(h - 1) 

We can see that these integrals have the same order behaviour as the corresponding 
homogeneous terms, that is, the Ki are 0(1) as t — > 0, for i = 1, . . . , 4 and k§ is 0(t~ c ) 
as t — > 0. 

Now since c > 0, (1551) shows that solutions to ( 1541) blow up as i 0. We now 
examine this divergence. 



5.2. Blow-up of the L q -norm 

We begin this analysis by determining a way to distinguish between those initial data 
which lead to diverging solutions to ( 1541) and those which do not. If we include both 
the homogeneous and the inhomogeneous parts, then we can label the five solutions to 
(1541 arising from Theorem 15.11 as 

0x0) = (l + Kl ,0,0,0,0) T , 
4> 2 (z) = (0,1 + K2 ,0,0,0) T , 
&(z) = (0,0,1 + /£ 3 ,0,0) T , 
fa(z) = (0,0,0,l + /t 4 ,0) T , 

$ s {z) = (o,o,o,o,r c + /£ 5 ) T , 

where the ^1,2,3,4 are finite as t — > and 05 is divergent. Given that ( 1541) has coefficients 
which are analytic on (0,ti], it follows that 0i_5 are analytic on (0,ii]. Thus these 
solutions provide a basis for solutions of ( 1541) on (0,£i]. Hence given any solution u(t) 
of (j5"4j) . there exist constants c?j, i = 1, . . . , 5 such that 

5 

u(t) = J2^i(t)- 

i=l 

Let S = L X (R, R 5 ). We consider solutions of (T441) with initial data in S. Given G 5, 
define Mo = f^u^ dp. We can define di(w-°>) via the existence of unique constants 
di, i = 1, . . . , 5 for which 

5 

M0 = 2Jrfi0(*l). 

i=l 

Define 

5' = G 5 : d 5 (u (0) ) = 0}. 
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Figure 3. The spread of the support of u. We illustrate the spread of the compact 
support of u from the initial data surface t\ to the Cauchy horizon. The growth of the 
support is bounded by in- and outgoing null rays starting from the initial data surface. 



This set corresponds one-to-one with initial data for (j44j) which give rise to solutions 
for which the corresponding solutions of (j54p are finite as t — > 0. 
We define the set complement of S' in S as S" = S — S'. 

Lemma 5.2 Given a choice of initial data u^ ' G S" , the solution u corresponding to 
this data displays a blow-up of its L x -norm, that is, 

lim I \u\ |i = oo. 

Proof We define u = J R udp where u is the solution of (jSJ) corresponding to E S" . 
It is immediately clear from fl55l) that 

lim \u\ = lim \u\ = oo. 
Then using the definition of the L 1 -norm of u, ||w||i, it follows that 



u\ |i = / \u\ dp > 



udp 



\u\, (56) 



which implies that \\u\\i — > oo as t — > 0. □ 

This divergent behaviour in the L 1 -norm of u could be attributed to the divergence 
of the support of u as it spreads from the initial surface t±, rather than to divergent 
behaviour in u itself (see Figure |3]for an illustration of this). We must therefore consider 
the behaviour of the spread of the support of u from the initial surface to the Cauchy 
horizon. In analysing the growth of the support of u, it is convenient to briefly return 
to the self-similar coordinate z. 

Lemma 5.3 Let the support of u be defined as vol[u](z) =p + (z) —p~(z) where 
p + (z) = sup{p : u(z,p) 0}, 
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p (z) = ini{p : u(z,p) ^ 0}. 

Then neglecting terms which remain finite on the Cauchy horizon, the support ofu grows 
as 

vol[u\(z) ~ — In \t\, t — > 0. 

Proof Define vol [it] (z) = p + {z) — p~(z) where 
p + (z) = sup{p : u(z,p) 7^ 0}, 

pGK 

p~(z) = inf{p : u(z,p) ^ 0}. 

pGM 

The support of u at some time z, supp[tt](z), will obey 

supp[u}(z) C [p~(z),p + (z)}. 
We define the L 9 -norm as usual, 







1/9 


rp+(z) 


\\u\\q 


1 \u(z,p)\ q dp 




I \u(z,p)\ q dp 








Jp-(z) 



1/9 



for 1 < q < oo. u has initially compact support which implies that vol[{t](^ ) = 
p + (z ) — p~(zo) < oo, where zo is the initial data surface. This initial support must 
grow in a causal manner; that is, the growth of p ± (z) must be bounded by the in- and 
outgoing null directions. From the metric ([T]), the in- and outgoing null directions are 
described by the relation dt/dr = ±e Uj/2 , which in (z,p) coordinates becomes 

dp [ h 
which results in 

dz{z±e v ' 2 )-\ (57) 

In handling this integral, we first substitute e u ^ 2 = S(z) — zS(z) (see (jSJ)), and then 
make the coordinate change y = S 1 ^ 2 (z). The resulting integral can be performed and 
results in 

4 

p+(z) = p + (z ) + 3 Uivt) ^ 1(1 + az) 1 ' 3 - y+\, 

i=l 
4 

p- (z) = p~ (z ) + 3 f- (%) In I (1 + az) 1/3 - y- | , 



where 
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and yf is the i th root of ±2a — 3k ± ak 3 + 3k 4 = 0. Therefore, the volume of the support 
of u grows as 

vo\[u](z) = p + (z ) - p-(zb) + (58) 

4 

3 £ (/+(%+) In 1(1 + azf' 3 - y+\ - /_(%) In |(1 + a^) 1 / 3 - % -|) . 
i=i 

Now, since in fl57|) there is a divergence at the Cauchy horizon when Zc = -e^ 2 (,2 c ), 
the above result must contain a Cauchy horizon divergence. Using the coordinate 
transformation y = S(z) 1 ^ 2 it is possible to show that in terms of y, the Cauchy 
horizon is at that y for which (y 3 — l)(3y — 2a) + 3ay 3 = 0, which is precisely where 
(1 + az) 1 ! 3 = —z(l + |az). When this holds, the first log term given above diverges. So 
in ( 158]) we have one finite term describing the initially finite volume, a second term which 
diverges on the Cauchy horizon and a third term which is finite everywhere. So overall, 
if we ignore terms which remain finite as t — > 0, then we can describe the behaviour of 
the volume as 

vo\[u}(z) In t->0. 

□ 

We will next need the L 9 -embedding theorem [16J, which we state as follows: 

Theorem 5.4 L q Embedding Theorem: Suppose Q C M n satisfies vol(Q) = 
J Q 1 dx < oo. For 1 < p < q < oo, if u E L q (Q), then u G L P (Q) and 

IMIp < [uoZ(fi)]p~9 \\u\\ q - (59) 

We are now in a position to show that the L 9 -norm of u diverges. 

Theorem 5.5 Given a choice of initial data E S" , the solution u of ( f^j ) 

corresponding to this data displays a blow-up of its L q -norm for 1 < q < oo, that 
is, 

lim \ \u\\ a = oo. 
Proof We set p = 1 in ( |59|) . This produces 

I Nli ^ [voiN(-2)] 1- ' I N !<?■ 

We know from fl56l) that ||"u||i > \u{z)\ ~ t~ c so 

r c 

1 — 1 

(volN(,z)) 8 

Now lim^o ^ c (l n (^)) 1 ^ = 0, since c > 0. Therefore, we can conclude that 
lim I NIL = oo. 

So the L 9 -norm of the solutions with initial data in S" blows up as t — > 0. □ 



< I N|„ (60) 
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We note that this behaviour will also hold for the four dimensional free evolution system, 
as the constraint is propagated (see Lemma [4. ID . We next prove two theorems which 
together show that this divergent behaviour is generic with respect to the initial data. 
In particular, we can show that the set of initial data which corresponds to solutions 
with divergent behaviour, S" , is open and dense in the set of all initial data, and that 
it has codimension 1 in S. 

Theorem 5.6 The quotient space of S' in S, S = S/S' , has codimension one in S. 

Proof A quotient space S = S/S' has dimension n if and only if there exist n vectors 
X(i), • • • , X(n) linearly independent relative to S' such that for every Ig5, there exist 
unique numbers C\, . . . , c n and a unique X' G S' such that X = Yli=i c iX% + X' [17J. So 
let X G S and let Xm be any element of S. To prove the result, we show that there is 
a unique value of a for which X^ a ' = X — aXp) G S'. Integrating over the real line and 
exploiting earlier notation, we have 

x^ = x — ax (i) 

5 

i=i 

We have X^ G S' if and only if d 5 (X) -ad 5 (X {1) ) = 0. Since d 5 (X {1) ) ^ - as X (l) G S 
- this occurs for a unique value of a. 

□ 

Theorem 5.7 S" is dense and open in S in the topology induced by the L l -norm. 

Proof To show that S" is dense in S, we must show that for any Ig5 and any e > 0, 
there exists some X" E S" such that the L 1 distance between X and X" is less than e, 
that is, 

WX-X"^ < e. (61) 

First, suppose that I G 5". Then (IBTl) is trivially satisfied by taking X" = X e S" . 
We therefore assume that X e S'. Now consider some ip(p) G C£°(R, R), such that 
ip(jp) > and J ip(p) dp = 1. We then set 

*" = * + >>mAv 

2 IWi)| 

where | • | indicates the maximum vector norm in R 5 , that is |</>5(£i)| = maxj|(0 5 (t 1 ))j| 
. Then 



\X-X"\\r = j ifo) 



\Mti)\ 

So we can explicitly construct the X" required to satisfy ( IBTl) . and thus S" is dense in 
S. 

| We note that if ^s(ti) = 0, we can replace 4> 5 with 5 = 4> 5 +J2t=i c i^>i^ f° r some constants Cj, chosen 
to guarantee that 4> 5 (ti) ^ 0. This does not affect the definition of S' or S" . 
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To show that S" is open in S, we must show that for all X" 6 S" , there exists an 
e > such that Bl(X") C S", where B\(X") indicates a ball of radius e in the Z^-norm 
centred at X" . We fix X" 6 5" and let X G ^(X") C 5. Then 

\\X-X"\\ l = I \X -X"\dp<e. (62) 
There exists unique constants Cj and (for z = 1, . . . , 5) such that 

x" = d\4>x + . . . + C?505, 
X = Cx4>x + . . . + C 5 5 . 

It follows from (162]) that |cj — di\ < c^e, for some a, depending on 1 _ 5 (t 1 ). Then by 
making e arbitrarily small, we can make the dj arbitrarily close to the q. We know that 
(I5 / since X" e 5"; therefore C5 7^ which implies that X & S". 

□ 

These two theorems, coupled with Theorem 15.51 suffice to show that the averaged form 
of the solution fl53|) displays a generic divergence of its L 9 -norm, where the term generic 
refers to the open, dense subset of the initial data which lead to this divergence. 

Remark 5.1 We note that if we define x := t c u and let x := j R xdp, then by 
multiplying (1331) by a factor of t c and taking the limit t — > 0, we find 

limx^O. (63) 
t-»o 

This will be used in the proof of Theorem 16. 101 The results of Theorems 15.61 and 15 . 71 tell 
us that the set of initial data which gives rise to (1631) is open and dense in L 1 (M, M 5 ). 

We note that from Theorem 15.51 alone we cannot conclude that the perturbation 
itself diverges on the Cauchy horizon. The reason for this is that one can easily imagine 
a function which has finite pointwise behaviour, but a diverging L 9 -norm arising from 
the spatial integration in (1331) . For example, a constant function is clearly pointwise 
finite, but has a diverging L^-norm. In the next section, we will determine the pointwise 
behaviour of the perturbation as the Cauchy horizon is approached and show that it 
diverges with a characteristic power of t~ c . 



6. Pointwise Behaviour at the Cauchy horizon 

So far we have determined the behaviour of the averaged perturbation u. In this section, 
we aim to show that the vector u has behaviour similar to that of u, that is, 0(1) 
behaviour in the first four components, and 0(t~ c ) behaviour in the last component. In 
this section, we will use Theorem 14.21 to provide us with smooth, compactly supported 
solutions u to ( 1441) . with a choice of initial data u(ti,p) G C^°(M, M 5 ). 

We begin by returning to the five dimensional symmetric hyperbolic system 

ni 1 ni 1 — * 

t^+A(t)^ + C(t)u = Z(t,p). (64) 
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Our strategy is to work with a scaled form of u, namely x := t c u. We can write an 
equation for x by using (164"1) . We find that 



Before presenting the results which determine the behaviour of x at the Cauchy horizon, 
we present a summary of various steps involved. 

• We begin by showing that x has a bounded energy throughout its evolution, 
including on the Cauchy horizon. Initially, we introduce the first energy norm, 
Ei[x](t), which is simply the L 2 -norm of x. In Theorem 16.14 we show that this 
norm is bounded by a term which diverges as the Cauchy horizon is approached. 

• We introduce a second energy norm, E 2 [x] (t) and in Theorem 16.21 show that it is 
bounded for t G [0, t*], for some t* sufficiently close to the Cauchy horizon. By 
combining theorems 16.11 and I6.2[ we can show that x has a bounded energy up to 
the Cauchy horizon; see Theorem 16.31 

• We use this to show that x itself is bounded in Corollary 16.41 However, there is 
no guarantee that x does not vanish on the Cauchy horizon. If this were to occur, 
then we would not be able to deduce any information about the behaviour of u 
from that of x. 

• We want to show that x is generically non-zero on the Cauchy horizon, for a set of 
non-zero measure. We can easily show that x := j R xdp is non-zero at the Cauchy 
horizon (see Remark 15. ip . If we could commute the limit t — > with the integral, 
then we could show that J R x(0,p) dp ^ 0, which would be sufficient, since then 
x over at least some interval on the Cauchy horizon. 

• In order to show that we can commute the limit with the integral, we turn to the 
Lebesgue dominated convergence theorem, which provides conditions under which 
one may do this. In order to meet these conditions, we must strengthen the bound 
on x (Lemma l6.5p . construct a Cauchy sequence of x values in L 1 (Lemmas 16.61 and 
16. 7p and finally apply the dominated convergence theorem. 

• So overall, we can show that x(0,p) ^ over at least some interval p G (a, b) on 
the Cauchy horizon. 

We begin with our first energy norm for x. 
6.1. Energy Bounds for x 

Since we expect u to diverge as t~ c , if we define x = t c u, then we expect that x should 
have a bounded energy in the approach to the Cauchy horizon. 

Theorem 6.1 Let u{t,p) be a solution of $J4\ ) subject to the hypotheses of Theorem \4-S\ 
Then there exists t* > such that for all t G (0, f ], the energy norm 



fir ( iv —* 



(65) 




(66) 
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obeys the bound 

Ei(t) < J, (67) 

for a positive constant fi which depends only on the initial data and the background 
geometry. 



[ 2t 2c - 1 {cu -u-u- C{t)u) - 2t 2c ~ l u ■ ^ + 2t 2c u ■ S dp. (68) 
Jr t op 



Proof We define E\{t) as in (166]) and take its derivative. We substitute f )65|) to find 

dEi _ f n4 2c-u^ n&^i.t 

Integrating by parts shows that the term containing duj dp vanishes due to the compact 
support of u and the fact that A(t) is symmetric. This leaves 

f 2t 2c -\cu-u-u-C{t)u) + 2t 2c - 1 u.tdp 1 (69) 
Jr 

on the right hand side of (1681 . We focus on the first term in (169]) . and introduce the 
constant matrix S, which transforms Co into its Jordan canonical form and which is 
listed in Appendix A. We use this matrix to show that 

cu ■ u — u ■ C(t)u = S T u ■ (cl — C(t))v, 

where v = S^u and C(t) = S' _1 C(t)S', so that C is the Jordan canonical form of C . 
Recall that C = diag(0, 0, 0, 0, c). Now since S T = S T S S^ 1 , we find that 

cu-u-u- C(t)u = v T S T S ■ (cl - C(t))v. 

Now for any matrix A, A T A is positive definite. We now wish to show that 
(v, S T S(cl — C{t)v)) > 0, and that equality holds iff v = 0. By using the form of S and 
the matrix C(t — 0) we find that the matrix S T S(cl — C(t) has four positive eigenvalues 
and one zero eigenvalue at t = 0. This indicates that it is positive semi-definite, but could 
still vanish along the direction of one of the eigenvectors. Specifically, it is possible for 
only the fifth component of v, v$ to be non-zero and for this dot product to still vanish. 
However, by the continuity of the matrices here, it follows that S T S(cI — C(t) has four 
positive eignevalues for t G (0, t*) for some t* > 0. Therefore, if (v, S T S(cI — C(t)v)) = 
in this range, it follows that we must have V\ = v-i = v% = v± = 0. But if these 
components vanish, then the only way to have (v,S T S(cl — C(t)v)) = is to have 
v 5 = too. We therefore conclude that if v\ = v<i = v 3 = f 4 = 0, we must also have 
v 5 = 0. Therefore, (v, S T S(cI - C(t)v)} = iff v = 0. 
So overall, we conclude that 

(v,S T S{d-C{t)v)) > 0. (70) 

for t G (0, t*) for some positive t*, with equality holding iff v = 0. We note that if v = 0, 
then u = Sv = too. 

We now assume that t G [0, t*) and using this information about the first two terms 
of (168]) . we can conclude that 



dt 



> / 2t 2c - L u-Zdp. 



Cosmic Censorship for Self-Similar Spherical Dust Collapse 



29 



We now apply the Cauchy-Schwarz inequality to show that 



> - t 2c - L u-udp- / t 2c -^-Hdp= — (E 1 (t) - J(t)), 

where J(t) := f R t 2c Yl • S. Integrating this from some time £ up to an initial time t\ 
(where < t < t\ < t*) will produce 

< hExih) + J J(r) dr. (71) 

We now examine the term J(t). We first recall the form of the source vector = 
h(t)g(p), where h(t) is an analytic function of t. J(t) can therefore be written as 
J(t) = t 2c h ■ hG, where G := J R g{p) 2 dp G R. Then the integral appearing in (fTTj) is 

-ti 

J(r)dT = Gtl c h-h (h-t), 



u 

where we use the mean value theorem [T5] to put the t-dependent terms outside the 
integral and t* G [t,ti\. The right hand side above is clearly bounded as t — > 0, and we 
can therefore conclude that 

Ei(t) < j, 

for a positive constant /i which depends only on the initial data and the background 
geometry. □ 

We note that f|66|) is in fact the L 2 -norm of x. Since E\(t) is bounded, it follows that 
x(t,p) G L 2 (R, R 5 ) for t G (0,t*). We also emphasise that this bound holds only for 
t > and does not hold on the Cauchy horizon. We now introduce a second energy 
norm E 2 [x](t) whose bound will extend to the Cauchy horizon. 

Theorem 6.2 Let u be a solution of [44\ )> subject to Theorem \4-S\ Define 

E 2 [x](t) := / x -x + (t- l)x\dp. (72) 

Then there exists some t 2 > and fi > such that 

E 2 [x^{t)<E 2 [x]{t 2 )e^- t \ (73) 

Here // is a positive constant that depends only on the components of the background 
metric tensor. 

Proof We begin by noting that the definition (172]) is equivalent to 

E 2 [x\(t) := / xl + x 2 2 + xl + xl + tx 2 5 dp. 

The factor of t is intended to control the behaviour of x$ as the Cauchy horizon is 
approached. In what follows, we will denote the (5, 5) component of the matrix A(t) as 
05 (i), the (5,5) component of the matrix C(t) as cs(t), and the remaining components 
of the fifth row of the matrix C(t) will be denoted c^j(t), where j = 1, ... ,4. We also 
recall that the first four rows of C(t) have 0(t) behaviour as t — y 0, and the fifth row 
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has 0(1) behaviour as t — > 0. We begin by taking a t-derivative of ( I72p and substituting 
from fl65l) to find 

dE 2 f ( A(t)dx C(t)-cl zti c _ x 



m j, 2 n * E(i ' p) J + ^ (74) 

where ^i=i c 5 j(t)xj appears in the source term of the x 5 equation due to the fact that 
the C(t) matrix is not diagonal. Now consider the terms 

A{t)dx x 5 dx 5 

- 2x -^^- 2(t - 1} T a5(t V 

This can be simplified to 

- 2 I xiax(t)— + x 2 a 2 {t)— + x 3 a 3 {t)— + x 4 a 4 (t) — 1 - 2x 5 a 5 {t) — , (75) 

where we used the fact that aj(t) = tdi(t) where aj(t) is 0(1) for % = 1, . . . , 4. After we 
insert f!75|) into the integral in f!74|) . it will vanish after an integration by parts, due to 
the compact support of x. Returning to ( 1741) . we are left with 

d -a¥= I 2 *' (-^^*+^(t,*)) +4+ (76) 
2{t - 1)^ ^-(c 5 (t) - c)x 5 + t c S 5 - Y,c 5j (t)x 3 ^J dp. 
If we now consider the terms 

-2X • (C(t) ~ Cl) f - 2(t - l)f (C 5 (t) - C) - 2(t - 1) £ ' 2 j !1 -'-rr:, 

3=1 

we notice that they can be rewritten as 

_2f • (^(^ ~ cI ) f _ 2 (c 5 (t) _ c ) x 2 _ 2 c 5j (t) Xj x 5 , 

where C(t) is a matrix got by replacing the final row of C(t) — cl with a row of zeroes. 
In other words, we can write C(t) as 



C(t) 



D 





where D is a 4 x 5 0(t) matrix. Equally, we could write C(t) = tC(t), where C(t) is an 
0(1) matrix. The fifth row of C(t) contains only zeroes, that is, C§i = for i — 1, ... ,5. 

We insert this into (j76|) and note that 2cx ■ x/t is explicitly positive definite. We 
also use the Cauchy-Schwarz inequality, in the form 

2x ■ £ dp > — / x-x + S- S dp, 
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to find that 
dt 



> / -2f-C(t)£-t c - 1 (£-£ + f-z) + a£ + 



-x 2 5 (c 5 {t) - c) + t c ~\t - - a£) - 2 c 5j {t)xjx 5 dp. 

3=1 

Now we let / equal the integrand on the right hand side above. We introduce 
Ir = I + /iIe 2 , where \i > is a constant and Ie 2 indicates the integrand of ([72]) . 
We wish to show that I R > 0. 
We can write Ir as 

I R = x- (-2C(t) + fd - t c " 1 I)f + (-2(c 5 (t) - c) + 1 - (t - ly- 1 

4 

+(t - l)n)xl - 2 Y c 5j (t)xjX 5 - • £ - (* - l)^ 1 ^. 

Let t = and note that c > 1 so that at t — 0, £ c_1 = 0. Then simplifies to 

4 



t=0 



x- (-2C(t = 0) + + (1-/^)^5 - 2^c 5j (t = 0)x i x 5 . 

3=1 



We can simplify matters by writing Ir 



x ■ Hx, where H can be written as 



t=o 



/ n \ 

// 
H= 00/i00 + K, 

^0 
\ / 

where if is a constant matrix, independent of //, which depends on the components of 
the matrix C(t) evaluated at t = and K 55 = 1. For Jr(0) > 0, we need H to be 
positive definite. This implies that all of the principal subdeterminants of H must be 
non-negative. So we require 



H + K n >0 






fi + K u 
K 21 


K 12 
^ + K 22 


>0, 




li + K u 


K 12 


K 13 




K 2 i 


+ K 22 


K 23 


>0, 


K 31 


K 32 


H + K 33 




H + Ku 


K 12 


K 13 


K 14 


K 21 


A* + K 22 


K 23 


K 2i 


K 31 


K 32 


H + K 33 


K u 




K i2 


K A3 


fi + K Ai 



>0, 
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\H\ 



ji + K u 


K 12 


K l3 


K u 


K 15 


K 2 i 


A* + K 22 


K 23 


K 2 a 


K 25 


K 31 


K 32 


H + K 33 


K M 


K 35 


K a 


K 42 


K 43 


li + Ku 


K 45 


K 51 


K 52 


K 53 


^54 


1 



> 0. 



These conditions produce a linear equation, a quadratic with leading /j 2 , a cubic with 
leading fi 3 and two quartics with leading //, all of which must be positive. We therefore 
pick a fi which is large enough to satisfy each of these conditions. 

We therefore conclude that at t — 0, Ir > 0. Now by continuity of the coefficients 
of x in Ir, it follows that there exists some t* such that Ir > in the range t G [0, t*]. 
We may therefore state that 
dE 2 



dt 



> -liE 2 , 



(77) 



in this range. We now integrate (177)) starting from some initial data surface t 2 € (0, t*) 
This results in 



M*2-t) 



E 2 [x](t) < E 2 [x]{t 2 )e 
which provides the desired bound for E 2 [x](t). 



(78) 
□ 



We note that the definition of E 2 [x](t), ( 1721) . is a sum of the L 2 -norms of Xi for 
% = 1,...,4 and the L 2 -norm of t 1//2 xs. Since we can bound E 2 [x]{t) in i G [0, t*], 
it follows that in this range, x { E L 2 (R, E) for i = 1, . . . , 4 and t 1/2 x 5 E L 2 (R, E). 

We next combine this result with Theorem 16. II to provide a bound on x which holds 
for the entire range of t. 



Theorem 6.3 Let u be a solution of $44\ ), subject to Theorem 
bounded by an a priori bound for t E [0, t\], that is, 

E 2 [x](t) < vEifflih), 

for a positive constant v. 



J3± Then E 2 [x}(t) is 
(79) 



Proof To prove this, we note that by definition, E 2 [x\{t) = Ei[x](t) + (t — l)Ei[x5](t). 
Therefore, using the bound on Ei[x\ from Theorem 15. II produces E 2 [x]{t 2 ) < Ei[x](ti). 
Inserting this into ( 173)) produces ( 1791 . where v = e^(t 2 — t). □ 



Corollary 6.4 Let u be a solution of ^J$ , subject to Theorem 
uniformly bounded in the range t E (0,t*). That is 



4-2. Then x = t c u is 



for i 



where the (3j, j 
the initial data. 



\Xi\ < (3i 

4 and 

\t l/2 x,\ < (3 5 

1, . . . ,5, are constants depending on the background geometry and on 
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Proof We first note that one of the effects of self-similarity has been to produce a 
differential operator on the left hand side of ( 165]) which has only ^-dependent coefficients. 
This means that the spatial derivative, x, p obeys the same differential equation as x, but 
with a modified source term. It follows that if we define E^x^it) := J t 2c u, p -u, p dp, 
we can bound this energy in an exactly similar manner to Theorem 15.11 Similarly, we 
can bound the energy E 2 [x^ p ] (t) using the same argument as that of Theorem 16. 21 
Now recall Sobolev's inequality which states that 

\v\ 2 < -x I W\ 2 + \v, P \ 2 dp, 
1 Jr 

for v G Co°(R, M 5 ). Applying this to x and using the bound (1791) (with the corresponding 
bound for x, p ) produces 

< Pi 

for i = 1 , . . . , 4 and 

\t 1/2 x b \ < (3 5 

where the (3j, j = 1, . . . , 5, are constants depending on the background geometry and 
on the initial data. □ 



6.2. Behaviour of x at the Cauchy Horizon 

Having established a bound on x through the use of energy norms, we now wish to 
determine the behaviour of x as t — > 0, that is, the behaviour on the Cauchy horizon. 
In particular, we must establish that x ^ there. To do so, we must first strengthen 
the bound on x by using the method of characteristics. 

Lemma 6.5 Let u be a solution of ftJ4\ ), subject to Theorem 4-2 Then x obeys the 
bounds 

\xi(t,p)\< lt t 1/2 , \x 5 (t,p)\< l5 , (80) 

for i — 1, ... ,4 and constants jj, j = 1, . . . , 5 which depend only on the initial data and 
the background geometry of the spacetime. 

Proof We begin by considering the first four rows of ( 1651) . which we write as 

dx ■ dx ' 

where Si(t,p) = t c Sj — Y^j=i j^i c ij{t) x j- Here a«(t) and c«(t) represent each entry 
on the main diagonal of the matrices A(t) and C(t) respectively (and we note that 
ai(t) = a 2 (t) = 0). Since the matrix C(t) is not diagonal, the term Y^j=i j^i c ij{^) x j 
represents all the off-diagonal terms, which we put into the source. 

We solve equations ( IHTI) along characteristics. The characteristics are given by 
p = pi(t) where 

dpi aiit) 



dt t 



ai(t)=>pi(t)=ni(t)+rH, (82) 
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where we use the fact that aj(£) is 0(£) as t — > 0, so that a,j(£) = £5j(£), where 5j(£) is 
0(1) as £ — > 0. 7Tj(t) = — J"* 1 di(r)dT and ^ = Pi(£i). On characteristics, flHTj) becomes 

t-^(t, Pi (t)) + (cj - c)^,^)) = Si(t, Pi (t)). (83) 
The integrating factors for these equations are given by e^ (£) where 

Jt r 

and if we Taylor expand the term inside the integral about £ = 0, we will find that 
e - ^ = (t/ti) c a>i(t), where a«(£) is an 0(1) term containing all terms other than the zero 
order term from the Taylor expansion. The solution to (1831) is 

xtfapi®) = (Pj C ^irs4 0) (Vi)-t C ^(t) T L_lV(r,^(r))dr, (84) 
\tij aitti) J t ai{r) 

where x\ is the initial data at £ = t\. Now we fix £ G [0, £i] and let pi G IBL Then using 
(1821) we can write dHl) as 



£ ^ aii(t) jo) 



rh T -c-i 

-t C ai(t) I -—Si(T,Pi(t) + 7Ti(r) - TTi(t))dT, 

Jt ®i{T) 

Now taking the absolute value of fl85|) will produce two integral terms (coming from the 
two terms in the source Si), which we label In and I^- In is given by 

hi = t c N(£)| / 1 ^—T c+1 \h t (T)\\g( Pi (t) + 7Tj(t) - ^{t))\dr, (86) 

Jt \<Xi{T)\ 

where we use the fact that Si(t,p) = £ c £j — Y^j=i j^i c ij{t) x j an d that Sj = thi{t)g{p) 
where the hi(t) terms are 0(1) functions as £ — > 0. We use the mean value theorem to 
evaluate this, and conclude that 

In = t c ^^Mt*)\\g(p t (t) + vr t (£*) - 7r*(t))l(*i - *), • (87) 

for £* G [£i,£]. That is, 

hi <^£ c + 0(£ c+1 ), (88) 

where ^ = sup te[0ttl] |/k(£*)||#fe(£) + n{t*) - ^{t))\t x . 
The second integral from ( 185]) is given by 



t\ T ~c-l 



1 



£>,(£)!/ j— ^ Ic^H^^p^+Tr^r) -<£))! ) (89) 



r -c-l 



+ F^TF ( I (r) 1 1 x 5 (r, Pl (£) + vr, (r ) - m (£) ) | ) dr. 
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To handle this integral, we first note that Cij(t) is 0(t), so that Cy(t) = tdij(t), where 
Cij{t) is 0(1) as t — > 0. We then use the bounds on x coming from Corollary 16.41 (as 
well as using the mean value theorem as before). This produces 

i -t° Mt)l ( V \c mi5 ( * r+1 - rc+1) ■ ic (t*)\3 ^ c+1,2 -t- c+112 ) 

\i=i jV* 
for £* e [*!,*]. That is, 

/o < ^ 1/2 + 0(t) (90) 

where u t = sup te[0itl] ^^|c i5 (t*)|/3 5 (-c + l/2)~ l . Combining ®, and (JHSD 

produces 

<Ti* 1/2 , (91) 

where the ji factors are constants depending on the initial data and the background 
geometry. 

Using (J9T]1 . it is also possible to improve our previous bound on x 5 (t,p). The 
equation which x§ obeys is 

t-^ + a 5 (t)-^ + (c 5 -c)x 5 = S 5 {t,p), (92) 
where S^{t,p) = t c £ 5 — ^- =1 c^j{t)xj. The characteristics for this equation are given by 



dp 5 a 5 (t) 



Ps(t) = n 5 (t) +i] 5 , 



dt t 

where 7r 5 (t) = — J^ 1 ^^-dr and 775 = Ps(ti). On characteristics, (192]) becomes 
dx ^ 

t-^-(t,p 5 (t)) + (c 5 - c)x 5 (t,p 5 (t)) = S 5 (t,p 5 (t)). (93) 
The integrating factor for this equation is given by e^ 5 (t) where 

Jt T 

and since c 5 (t = 0) = c, we see that e~^ 5 (t) is an 0(1) function as t — > 0. The solution 
to (El is 



x 5 (t,p 5 {t)) = e-&4 0) (P5(*) - tt(0) (94) 

t /•*! e &(V) / A A 

- e_?5 y — ^ ^ c S 5 (r,p 5 (t) + 7r 5 (r) - 7r 5 (*)) - 22c 5j {T) Xj {T,p 5 {t) + 7r 5 (r) - 7r 5 (t))J dr. 

The integral above contains two terms, 

h = e~ €s J'' ^^r c S 5 (r,p 5 (t) + 7r 5 (r) - 7T 5 (t))dr, (95) 

and 

j 2 = e -€ 5 / I Vc^rj^fr.ftff) + 7r 5 (r) - 7r 5 (t)) ] dr. (96) 
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We recall that E 5 (i,p) = hs(t)g(p) and use the mean value theorem to show that 

Ji = e-^(t*)h 5 (t*)g(p 5 (t) + 7r 5 (r) - n s (t))^Q-, (97) 

c 

for t* G [ti,t\. For the second integral, we use the mean value theorem and the bound 
( 19T|) which arises from the first part of this theorem. This produces 

h = \2e-^eHt*) p^ 5j (t*)j/J (t 1 / 2 - (98) 

for t* e [tx,t]. Combining §7$, ([98]) and flM} produces 
\x 5 (t,p)\ < 75, 

where 75 is a constant depending only on the initial data and the background geometry. 

□ 

The next lemma shows that we can bound t 1//2 Xj,£. We use this in Lemma ISTTl to construct 
a Cauchy sequence of x- values in L 1 . 

Lemma 6.6 Let u be a solution to [44\ )> subject to Theorem \4-S\ Define x '■= dx/dt. 
Then 

\t 1/2 Xi(t,p)\ < V i, (99) 

for i — 1, ... ,4, where r\\ is a constant depending only on the background geometry and 
the initial data. 

Proof Define \ := dx/dt. By differentiating ( 165]) . we see that x obeys 

*lt + m % + {C{t) + ° m = * ( '' p) ' (1 ° 0) 
where a = (t c S), 4 —A, t §| — C, t x. For i = 1, . . . , 4 this becomes 

t^ + a i (t)^ + (c i (t) + l-c) X i = cr i (t,p)- £ Cij (t) Xj = Si(t,p), (101) 

where cij(t) and Cj(t) label the diagonal elements of the A(t) and C(t) matrices appearing 
in (llOOp and since C(t) is not diagonal, the %(£) label the off-diagonal elements which 
we put in the source term Si(t,p). We note that ai(t) = ct 2 (t) = 0. 

As in Lemma 16.5] we solve (llOip on characteristics (and these are the same 
characteristics which appeared in Lemma 16 .5 j) . By a similar method to that which 
lead to ( |85l . we find that the solutions to (110 ip can be written as 

Xi (t,p) = a % {t) (P\ X f\pi ~ 7Ti(t)) - (102) 



ai(t)t c 1 / — —Si(T,pi - -Kiij) + TTi(t))dr, 



where xf\pi ~ ^(O) is the initial data and the integrating factor is e~^(t) = 
t c ~ 1 t\~ c ai(t). As in Lemma 16.51 we Taylor expand the integrating factors and put 
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all higher order terms into 0(1) functions aa(t). We note that c — 1 > 0. Taking an 
absolute value of f l 1 2 j) produces two integrals (coming from the two terms in Si(t,p)), 
which we label 

I a {t,p) = \ai(t)\t c - 1 [ 1 -l—la^p^r) + tt^t) - Ki(t))\dT, (103) 



Ht,p) = \a l {t)\t c - 1 —— TCiA T )xA T MT)+n{T)-Ki{t))dT, (104) 



.. , Oij\T) 

and 

where we recall that Cy(t) is 0(t), so that Qj(t) = tcij(t), where Cjj(t) is 0(1). Now to 
deal with J il( we note that Gi(t,p) = (t c £j), t —a ilt ^ — C, t x\%, where C, t x\i indicates the 
% th row of the matrix C, t x. Using the bounds flHDl on x (and note that C, t x\i includes 
x&), we can see that overall o~i(t,p) is 0(1) in t. It is also C°° in t, so we can apply the 
mean value theorem to find 

a 8 (r)a,(r,^(r) + - m(t)) (t- c+l - r c+1 ) 



Ia(t,p) = t 



c-1 



an(t*) (-c + 1) 

for t* G [t\,t\. We can abbreviate this by writing 

in{t,p) <m, (105) 

where rji is a constant depending on the initial data and the background geometry 
(it inherits this dependence from the bounds on x entering into o~i). For ija, we note 
that from ( l65l) . we can deduce that \t\\ is 0(1) (actually, the only 0(1) term is the 
a^{t)dx^/ 'dp term) and it is also C°° in t, so as before, we apply the mean value theorem 
to find 



5 



for £* G [ti , t] , which we can summarise as 

Ii2(t,p)< m t. (106) 
Combining f |102j) . ( I105|) and f ll06|) produces 

k<,t(*,p)| = |Xi(*,P)l<^, (107) 
for i = 1, . . . , 4, where we neglect higher order terms. 

□ 

We next use this result to show that we can define a sequence x^ of if-values which is 
Cauchy in L 1 (R, R 5 ). 

Lemma 6.7 Lei {i*™- 1 } fre a sequence oft-values in (0, t\\ with lim^o^ — 0. For each 
n> 1, define x^(p) = x(t^ n \p). Then {x^} is a Cauchy sequence in L : (R, R 5 ). 
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Proof We define x^ := x(t^ n \p), where the sequence {t^}^Lo t en ds to zero as n — > oo. 
The mean value theorem produces 

\x{t (m \p) - x{&\p)\ = \x, t {t*)\\t (m ^ - tW\, (108) 

for some t* e (t (m) ,t (n) ). For % = 1,...,4, we can use the bound from Lemma 16.61 and 
integrate with respect to p to give 

\\x i (t ( - m \p)-x i (t^,p)\\ 1 < \tW-tW\rn / dp, 



pi 

where p\ = max peR (supp[£;(7;( m ),p)], supp[xj (t (n) ,p)]) and 

p 2 = min pgR (supp[a; i (t^ m \p)], supp[xi(t^,p)]). From Lemma |5.3[ we know that 
supp[xj(t,p)] satisfies supp[xj(£,p)] ~ ln(t) + /x, where \i is a constant which represents 
terms that remain finite as t — > 0. Now since < t (m) < i (ri) < t x , the largest support 
is that at t^ m \ so that 

1 1 Xi (t (m) ,p)-x i (&\p)\\ 1 < |t (m) - t (n) 1 77! (ln(t (m) ) + /i) . 

We can take the n — > oo limit above and see that x\ n \ for i = 1, ... ,4, is a Cauchy 
sequence with respect to the L 1 -norm. 

To show that Cauchy sequence in L 1 we use a different tactic. As in (194|) . 

we write the solution to the 25 equation of motion as 

x 5 (t,p) = e-^(t)xf\p - tt 5 (*)) - F(t,p), (109) 

where 

F(t,p)=e-« 5 (t) f 1 —S 5 (T,p + 7r 5 (T)-7r 5 (t))dT. (110) 

Here e^ 5 (t) is the integrating factor, and £5(7;) = — J. 1 C5 ^~ c dr, 7r 5 (t) = — J 1 ^^-dr 
and cs(i) and a^t) are the (5, 5)-components of the C(t) and A(t) matrices respectively. 
S 5 (t,p) is the source term, S 5 (t,p) = t c H^(t,p) — J2j=% c 5j(t) x j{t>P)> where c 5 j(t) are the 
components of the fifth row of the C(t) matrix appearing in fl65l) . Using (11091) we can 
write 

\x 5 (&\p) - x 5 (t (m) ,p)| < (111) 
\e^(t^)xf\p-7r 5 (t^)) - e ^(t^)xf\p-7r 5 (t^))\ + \F(t^ ,p) - F{t™ ,p)\. 

Now if we suppose that t {n) and 7j (m) are very close to the Cauchy horizon, t = 0, then 
tracing the characteristic back to the initial data surface, we can see that it will intersect 
the initial data surface outside the compact support of the solution. That is, there exists 
some N such that for n,m> N, xf\p — Tr 5 (t^)) = xf\p — ir^t^)) = 0. 

Now for the second term in ( 11 111) , we use the mean value theorem to show that for 

t* e [t( n \t^\, 

\F(t {n \p) - F{t {m \p)\ < ^-(t*,p) \t {n) - t {m) \. (112) 

at 

We can easily calculate 
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Now using the bounds (1801) on x (and the fact that S 5 = h 5 (t)g(p) where h 5 (t) is 0(1) 
as t — y 0) we can show that 

where fi is a constant depending only on the background geometry and the initial data. 
Combining this with ( I112p produces 

\F(&\p) - F(t (m) ' p )| < ^(t*)- 1/2 |t (n) - t (m) |. (114) 
Returning to fillip , assuming n,m > N and using (I114p produces 

\x 5 {t (n) ,p) - x 5 (t (m) ,p)| < /i(t*)" 1/2 |t (n) - t (m) |. 
Finally, we take the L 1 -norm to find 

\\x^ n \p) -x 5 (t (m) ,p)||i < fi(tT l/2 \t (n) -t (m) \ / ^, 

where p x = max pGR (supp[x 5 (t^ m \p)], supp[a;5 (t (n) ,p)]) and 

^2 = Hiin pe ]K(supp[a;5(t^ m '\p)], supp[xs(t^,p)]). From Lemma |5.3[ we know that 
supp[x 5 (t,p)] satisfies supp[xi(t,p)} ~ ln(t) + /i, where \i is a constant which represents 
terms that remain finite as t — > 0. Now since < t (m) < < ti, the largest support 
is that at t^ m \ so that 

\\x 5 (&\p)-x 5 (t^\p)\\ 1 < v(tT 1/2 \t {n) -t {m) \(Mt^ +f ,). 

and taking the limit n — > oo, we see that x^ is also a Cauchy sequence in L 1 . So we 
can conclude that x^ is a Cauchy sequence in L 1 (IR, IR 5 ). 

□ 

Now with Lemma [6.71 in place, and since we know that x G L 1 (IR, IR 5 ) for t G (0, ti], we 
can show that x does not vanish on the Cauchy horizon. To do this, we will make use 
of two theorems from real analysis, which we state here. The proofs of both theorems 
are standard; see [19] for details. 

Theorem 6.8 Let 1 < p < oo. Suppose Vt is some set, Q C R n , and let /W, % = 1, 2, . . . 
be a Cauchy sequence in L P (Q). Then there exists a unique function f G L p such that 
\\f - /||p as i -> oo ; that is, /W converges strongly in the L p -norm to f as 
i — > oo. 

Furthermore, there exists a subsequence /fr), /(«),..., h < i 2 < ... and a non- 
negative function F G L P (Q) such that 

• Domination: |/^ fc )(a;)| < F(x) for all k and for a dense subset of x EVt, 

• Pointwise Convergence: lim^oo p %k >(%) = f(x) for a dense subset of x G O. 

Theorem 6.9 (Lebesgue's Dominated Convergence Theorem) Let {/^} be a sequence 
of summable functions which converges to f pointwise almost everywhere. If there exists 
a summable Fix) such that \f^'{x)\ < F(x) Vz, then \f(x)\ < F(x) and 

lim / f®(x)dx= / f(x)dx, 
'-► ° Jn Jn 

that is, we can commute the taking of the limit with the integration. 
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Theorem 6.10 Let u be a solution to (J$, subject to Theorem\4.^\ Then x(t,p) = t c u 



does not vanish as t — > 0, for a generic choice of initial data. Here the term generic 
refers to the open dense subset of L 1 initial data which lead to this result. 

Proof Since fei 1 for t G (0, ti], the proof of this theorem follows by an application of 
the two theorems from analysis quoted above. Consider the sequence x 1 ^ := x(t^ n \p), 
where {t^}^ =0 tends to zero as)i-> oo. In Lemma [6.71 we showed that this sequence 
is Cauchy in L 1 . Theorem 16.81 therefore provides for the existence of a dominated 
subsequence of x^ n \ In particular, by applying this theorem we may conclude that 
there exists a non-negative H G L 1 (IR, R) and a unique h G L 1 (IR, R 5 ) such that 

\x M \<H(p) Vm, 

and — h\\\ — > as m - » oo. Next we apply the Lebesgue dominated convergence 

theorem (Theorem 16.91) to the dominated subsequence x^ nm ' . This produces 

lim / x(t { - rim \p)dp= / x(0,p)dp, (115) 
m ^°° Jr Jr 

where we identify h G L^R, IR 5 ) with x(0,p). Q 

Now if we recall Remark 15.11 which indicated that lim^o j^^dp 7^ 0, we can 
conclude that lim, n _ >00 j R x{t^ nm \p) dp ^ 0. Combining this with f 1 1 1 5 j) produces 

x{0,p)dp^0, (116) 

which implies that there exists an open subset (a, b) on the Cauchy horizon such that 
x{tiP) for p G (a, b). We note that (11161) holds generically since lim^o f R ^dp ^ 
for a generic (that is, open and dense in L l ) set of initial data (see Remark 15. II) . □ 

We conclude that x := t c u exists and is non-zero on the Cauchy horizon for p G (a, b), 
for a general choice of initial data. This in turn tells us that the perturbation u diverges 
in a pointwise manner at the Cauchy horizon, with a characteristic power given by t~ c . 

Theorem 6.11 There exists an open and dense subset of initial data G L 1 (R, R 5 ) 
such that the solution u of §6J$ corresponding to this initial data blows up as t — >■ on 
an open subset p G (a, b), that is 

\imu(t,p) = oo, \/p G (a, b). (117) 

Proof It follows immediately from Theorem 16.101 that u blows up as t — > on an 
open subset p G (a, b) for a choice of C^°(R, R 5 ) initial data. Recall Remark 15.11 which 
indicates that x = J R xdp ^ for a generic choice of initial data from L 1 (R, R 5 ). We 
can therefore extend the results of Theorem 16.101 to a choice of initial data taken from 
an open and dense subset of L 1 . We conclude that for such initial data 

\imu(t,p) = oo, Vp G (a, b). 

f That is, x(0,p) is defined on a dense subset of M. by the second result of Theorem 16.91 It suffices to 
take any bounded extension to "fill in" the definition of x(0,p) on the remaining set of zero measure. 
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7. Physical Interpretation of Variables 

So far, we have established the behaviour of u as it approaches the Cauchy horizon. We 
now wish to provide an interpretation of these results in terms of the perturbed Weyl 
scalars, which represent the gravitational radiation produced by the metric and matter 
perturbations. In this section, we will use the coordinate system (u,v,8,(f)), where u 
and v are the in- and outgoing null coordinates (see (TTTT) for their definitions), rather 
than (z,p, 8, <ft) coordinates. This is a useful choice of coordinate system to make when 
considering the perturbed Weyl scalars. We follow throughout the presentation of [21J. 

For the even parity perturbations, only two of the perturbed Weyl scalars, 5*o 
and 6*4, are identification and tetrad gauge invariant (see [21] and [22]). This means 
that if we make a change of null tetrad, or a change of our background coordinate 
system, we will find that these terms are invariant under such changes. We note that 
5*o and 6*4 represent transverse gravitational waves propagating radially inwards and 
outwards. These terms are given by 

5* = ^l A l B k ABl = ^- 2 n A n B k AB} 

where Q and Q* are angular coefficients depending on the other vectors in the null 
tetrad, and on the basis constructed from the spherical harmonics. The ingoing and 
outgoing null vectors l A and n A are given in f )12p . The term 

6P-! = |<W^ 4 | 1/2 (118) 
is also invariant under spin-boosts, and therefore has a physically meaningful magnitude 



Theorem 7.1 The perturbed Weyl scalars <5\J/ an d 8^4, as well as the scalar 5P_i, 
diverge on the Cauchy horizon. 



Proof We begin by writing the tensor k AB in (w, v) coordinates as 
k AB = 



T](u,v) u(u,v) 

u(v,v) X(v,v) 



The condition that k AB be tracefree results in u(u,v) = 0. In (u,v) coordinates, the 
perturbed Weyl scalars become 

Q Q* 
5 * 0= 2^p2 V(u,v), 8Vt = —\(u,v), 

where the factor of B(u, v) is the same factor which appears in (TI2"]) . Now by performing 
a coordinate transformation, we can write a(z,p) and (3(z,p) (the components of k AB 
in (z,p) coordinates) in terms of rj(u,v) and X(u, v) and by so doing, we can find 5*o 
and 6*4 in terms of a(z,p) and f3(z,p). We find 

5* = F(z,p)(a(z,p) - f:\z)f3(z,p)), (119) 
6*4 = G(z,p)(f + (z)a(z,p) - P(z,p)), 
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where the coefficients F and G are given by 

Q fl ( f- 
f- - A 



2r 2 v 2 \f + -f_ 



,-2 



(recall that B(u,v) 



F ^ Z ' P) 2r 2 B 2 u 2 

and we note that F ~ r -4 and G 

Now, if we retrace our steps through the first order reduction in Section HI we find 
that us(z,p) contributes to a(z,p), fl(z,p), k(z,p) and its first derivatives. In particular, 
the pointwise divergence of u$(z,p) on the Cauchy horizon produces a similar divergence 
in these terms. We can write a and fl as 
S 



a(z,p) 



w 4 (-l + S)+ u 5 (l + S) + 2Ui(l + 5) 



(120) 



1-S 2 

/3(z,p) = (u 4 (-l + - 5 + zS) + 2z(l + S) Ul + u 5 (z -zS + S(l + S)) 

1 + S 2 v 

where S(z) = (1 + az) 2//3 is the radial function. Combining (j!19p and (j!20p produces 

5^o = F(z,p)(fl 1 (z)u 4 (z,p) + fl 2 (z)u 5 (z,p) + fl 3 (z)u 1 (z,p)), 
= G(z,p)(fl 4 (z)u 4 (z,p) + p 5 (z)u!(z,p) + fl e (z)u 5 (z,p)), 



where 



fli(z 

fll{z 

flz{z 

fli{z 

fl 5 (z 
fleiz 



2(-l + S)SS 2 
1-S 4 ' 
1 z-zS + S(l + S) 




(z - S + zS)(l + S 2 ) 

z(l + S) 
1-5 " (z - S + zS)(l + S 2 ) 

2 (~ 1 + S \ s + zS{-l + S)), 
1-S 4 

2S 2 (1 + S) 2 
1-S* ' 

' S {-S(l + S)(-2SS - z + S 2 z)). 



1-S 4 

So 5^o and 5^ 4 depend on u§ and therefore they diverge as the Cauchy horizon is 
approached. Similarly, b~P-\ diverges on the Cauchy horizon, as it depends on 5^o and 

To construct a gauge invariant interpretation for the matter term T(z,p), we note 
that by comparing the GS terms (|20|) and f l24|) to fl32|) . we find that 

T A = p(z,p) u A (T(z,p) +g(p)) 
which we contract with the background dust velocity u A to find 

u A T A = -p(z,p)(T(z,p) +g(p)) 
where u A T A is a gauge invariant scalar. 
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8. Conclusion 

We have considered here the behaviour of even parity perturbations of the self-similar 
LTB spacetime in the case where there is a naked singularity. In particular, we have 
examined their evolution in the approach to the Cauchy horizon associated with the 
naked singularity. We first identified the fundamental system of PDEs which govern 
the evolution of the even parity perturbations but chose to use the five dimensional 
system in preference, as its symmetric hyperbolicity is a very useful property. We next 
considered an averaged form of the perturbation variable and showed that it displayed a 
generic divergence in its L 9 -norm. The divergence occurs for all angular number I. Here 
the term generic refers to the fact that the divergence occurs when solutions evolve from 
a set of initial data which is open and dense in the set of all initial data. This result is 
not by itself sufficient to ensure that the perturbation diverges on the Cauchy horizon, 
as it is perfectly possible for a pointwise finite function to have a diverging L^-norm. 
However, this result gave us a useful ansatz for investigating the asymptotic behaviour 
of the solution itself in the approach to the Cauchy horizon. 

To determine the pointwise behaviour of the perturbation in the approach to the 
Cauchy horizon, we introduced a scaled version of the state vector, x := t c u. We used 
a series of energy methods which produced bounds on x and its p-derivatives in the 
approach to the Cauchy horizon. By combining these results with the formal solutions 
for x arising from the method of characteristics, we were able to establish that fei 1 . 
Finally, this allowed us to establish a theorem which demonstrated that x is non-zero 
on the Cauchy horizon over some interval. This result pertains for all angular numbers 
I. 

Iguchi, Harada and Nakao [23] studied the behaviour of the quadrupole mode (/ = 2) 
of the even parity perturbations of the LTB spacetime. They numerically solved the 
linearised Einstein equations and found that this perturbation diverged on the Cauchy 
horizon. In one sense, our results are a generalization of theirs, in that this method 
allows us to treat all perturbations. However, the extra symmetry of self-similarity is 
needed in order to apply our methods. 

We note a potential issue with this work. Our perturbations are at linear order. It 
follows that it is somewhat strange to conclude that they diverge on the Cauchy horizon, 
as they are therefore far too large to remain at linear order. However, this result still 
indicates that this spacetime is not stable to perturbation. 

In previous work, we demonstrated that the odd parity perturbations of this 
background spacetime are finite for all I, where finiteness was measured with respect to 
initial values of a natural energy norm for the odd parity system. Taken as a whole, 
this work supports the hypothesis of cosmic censorship, in that one should expect 
perturbations on a naked singularity spacetime to diverge as the Cauchy horizon is 
approached. 

The background spacetime investigated here is of course not a serious model of 
gravitational collapse, as at the very least, it ignores the effects of pressure during the 
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collapse. A natural extension of this work would therefore be to consider the self-similar 
perfect fluid model, which contains a naked singularity for a wide range of the equation 
of state parameter. The study of the behaviour of perturbations in this spacetime would 
be a very interesting application of the methods developed here. 
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Appendix A. Equations of Motion 



We list here the various matrix coefficients and source terms omitted in sections HI |5] and 
EJ We neglect to list those terms whose exact form is not important. We note in what 
follows that S = S(z) = (1 + az) 2 / 3 is the radial function, q = q(z) = \/A-k(—zS + S)S 2 
is the density function and ■ = ■§:■ 

The five dimensional system takes the form 

du ~ Ou ~ 

— + A(z)— + C(z)u = a(z,p) 

where the matrix coefficient A(z) is given in Section HI C(z) is given by 



C 



( Cn 


C12 


Cl3 


C14 


Cl5 \ 


C21 


C22 





C24 


C25 


C31 


C32 


C33 


C34 


C35 


C41 





C43 


C44 


C45 


V c 5 i 





c 53 


C54 


c 55 / 



where 



Cll 



-3z(-l + S)S 2 - S 2 S + S(3S 2 -zS + S{-3 + 2zS)) 



S{-1 + S){S - zS) 



-l + S 



■1 + S 



C12 = 


2S ' 
-1 + 5 


Cl3 = 


25 ' 
-1 + 5 


C14 = 


2S ' 
8nqS 


C15 — 


25 ' 

_4 


C 2 1 = 


-1+3' 


C22 — 








q 




AirqS 




AnqS 


C24 = 


1 + 5' 


C25 = 


-1 + 5' 





ni 

° 31 ~ z5(-l + 5)(5-25) 2 ' 

ni = 2(2^ 3 5 4 - 2^ 2 55 2 (35 + z5) + 5 3 (-25 + z 2 S (3) ) 
+zS 2 (QS 2 + z 2 S 2 + zS(2S - ^5 (3) ))), 
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~ 2*5(1 + S) ' 

n 2 = zS{2 -L-L 2 + 2S + 2S 2 ) - AS 2 S + S{-2 + L + L 2 -2S 2 - 2zS 
+S(-2 + AzS)), 



2zS(-l + S)' 

n 3 = zS{-2 + L + L 2 + 2S - 2S 2 ) + AS 2 S - S{-2 + L + L 2 - 2S 2 + 2zS 
+S(2 + AzS)), 

4{l + S){-zS 2 + S{S + zS)) 
° 41 ~ S(-l + S)(z-S + zS) ' 

(l + 3)(S-zS) 
S(z-S + zS) 

3 + {z + S)S + S(3 + zS) 

° 44 ~ {l + S)(z-S + zS) ' 

(l + S)(S-zS) 

C45 — = — 1 

S(z-S + zS) 

4{-zS 2 + S{S + zS)) 
° 51 ~ S(z + S-zS) ' 

j-l + S){S-zS) 
° 53 ~ S(z + S-zS) ' 

(-l + SVS-zS') 

C54 = : , 

S(z + S — zS) 
3+(-z + S)S + S(-3 + zS) 
° 55 ~ (-l + S)(z + S-zS) 




C33 = - + 







2S_ 

S' 




f(t)g(p) where /(t) 



8nq(l + S)(S - zS) 
S(z-S + zS) 
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The coefficients appearing in (|42p are given by 

AS{-6z 2 SS 3 + 2z 3 S A + AnzqS 2 (S - zS) 2 + S 3 {-2S + zS) + zS 2 {6S 2 - zSS + z 2 S 2 )) 



9i{z) = 

02(2) = 
93(z) = 

9a(z) = 

9b(z) = 
where 

7711(2) 



2S(2S 2 + 2z 2 S 2 + zS(-AS 

-2S 2 {S + z{-S + zS)), 
m\(z) 

"T+l' 

2S 2 (S-zS), 



(-1 + S)(-S + zS) 
z§)), 



95{z) 
9i{z) 



m 2 (z) 

32irqS(S - zS) 2 , 



m 2 (z) 



S(z 2 S(l + S)(-2 + L + L 2 -2S- 2S 2 ) + S 3 (87rzq - AS) 

-zS{-2 + L + L 2 -AS 3 - 2zS + 2S{-3 + L + L 2 + zS) + S 2 (-6 + AzS)) 
+S 2 {-2 + L + L 2 -2S 2 + AzS + S{-2 - 8nz 2 q + 8zS))), 
S(-z 2 (-l + S)S{-2 + L + L 2 + 2S- 2S 2 ) + AS 3 (2itzq + S) 

-zS{-2 + L + L 2 + AS 3 + 2zS - 2S{-3 + L + L 2 - zS) - 2S 2 {3 + 2zS)) 
-S 2 {-2 + L + L 2 - 2S 2 - AzS + S{2 + 8itz 2 q + 8zS))). 

The coefficient matrix E(z) appearing in (143jl is given by 

/ 1 \ 



„, 3(1 + a ? )1/3 ,. 





E(z) 



\ 



i+az— 

3^(l+az) 1 /3 





3+a^+3z(l+az) 1 / 3 / 



so we can explicitly see that this system is not symmetric hyperbolic. 

Finally, we present the similarity matrix S (which appears in Section 16.11 and is 
used in the proof of Theorem 16 .ip . a constant matrix which transforms the zero order 
term of the Taylor expansion of C(i), C(t — 0), into Cq. It is given by 

/ si s 2 s 2 \ 
1 



S 



10 
10 
\ 1 1 J 



where 



si 



S2 



S{z c ){3 + {-z c + S{z c ))S{z c ) + S{z c ){-3 + z c S(z c ))) 
4(-l + S{z c )){z c S{z c ) 2 - S{z c ){S{z c ) + zS(z c ))) 
(-l + S(z c ))(-S(z c ) + z c S(z c )) 
A(-z c S(z c ) 2 + S(z c )(S(z c ) + z c S(z c ))) ' 
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Appendix B. The Behaviour of c 



In the statement of Theorem 15.11 we looked that the Jordan canonical form of the zero 
order term in the Taylor expansion of the matrix C(t). We claimed that the only non- 
zero eigenvalue of this matrix was in the range c G (3, +oo), for a G (0,a*). Here we 
prove this claim. 

In the previous section, we presented the coefficients of the equation of motion for 
u. Recall that the coefficient of u was given as C(z), where in terms of the notation of 
(j44p . C(t) = tC(z). After we put C(z) in Jordan canonical form, we find that the only 
non-zero term is the (5, 5) entry, which is given by 

3 + (S-z)S + S(-3 + zS) 

= ; : . B.l 

(1 - S)(z + S - zS) V ; 

Now define h(z) := z + S — zS. On the Cauchy horizon, h(z c ) = by definition. We 
can use this to Taylor expand the numerator of ( IB. II) and to simplify the zero order 
term. We can also Taylor expand the factor of h(z) which appears in the denominator 
so that h(z) = h(z c )(z — z c ) + 0((z — z c ) 2 ) = (1 — z c S(z c )) + 0((z — z c ) 2 ). Inserting 
these expansions into (IB. 1[) produces 

5 55 = ( 2 + K—^l + 0(1), (B.2) 



1 - z c S(z c 



as z — > z c . So we must determine the value of h(z c ) = 1 — z c S(z c ) = 1 + |a z c (a)(l + 
az c (a))~ 4 / 3 , where we use ((7|) to prove the last equality and we write z c = z c (a) to 
emphasise that the location of the Cauchy horizon depends on the value of a, the 
nakedness parameter. 

We first analyse how z c depends on a, before using this information to determine the 
behaviour of h(z c ). We define X(z) = (l + az)h(z) 3 , so that X(z) = z 3 (l + az) + (l + ^z) 3 . 
The Cauchy horizon corresponds to the first negative root of the quartic equation 
X(z) = 0. This root must lie in the interval z G (—^,0), since z = — - corresponds 
to the singularity (see ©). We know that a root exists for a G (0, o*), where a = a* 
corresponds to a double root of A, where X(z c (a*)) = X(z c (a*)) = 0. Let z* := z c (a*). 
We can easily show (using A = 0) that 

(l + ~z)A = z 2 (3 + Aaz + \a 2 z 2 ), 



3 



A=0 3 



which implies that a*z* satisfies the quadratic 

(a*) 2 (z*) 2 + 12a*z* + 9 = 0. (B.3) 
By solving this for a*z* and picking the larger root, we have 

z* = l(-6 + 3>/3). 
a 

Next, we consider the dependence of z c on a. If we differentiate the condition A = 
with respect to a, multiply by (1 + az) and use A = to simplify, we find that 

dz 

(a 2 z 2 + Ylaz c + 9)— - = 2az 3 c . (B.4) 
da 
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We note that the coefficient of dz c /da vanishes at a = a* (see lB.3p and is positive in the 
interval a £ (0, a*) (we can see this by calculating its roots and noting that for the range 
of allowable a, we are always above the larger root). Now since z c < 0, this implies that 
dz c /da < for all a £ (0, a*) and 

hm — — = — oo. 

a^a* dd 

So to summarise, we know that z c (0) = —1 (see the definition of A), z c (a*) = z* and z c 
is monotonically decreasing from —1 down to z* as a increases from to a*. 

We now determine the range of h(z c ). Define u(a) = h(z c (a)). Then u(0) = 1 (since 
u(a) = h(z c ) = 1 + |a 2 z c (a)(l + a,2 c (a))~ 4//3 ) , and by definition u(a*) = (recall that 
A(a) = (1 + az)h(a) 3 ). A straightforward calculation (using (1B.4jl ) shows that 

du 4 . A,n (3 + 2az c ) 

da 9 3 + 4az c + ^a l z l c 

Since a > 0, z c < 0, 1 + az c > and 3 + 4az c + \a?z 1 c > on a £ (0, a*), it follows that 
du/da < for a £ (0, a*). So w(a) = h(z c (a)) £ [0,1] and w decreases monotonically 
from it = 1 at a = to u = at a = a*. It follows that £ (1, oo) for a £ (0, a*) with 

lim ——- = 1, lim — — = +oo. 

a^o+ u{a) a->(o*) _ tt(a) 

We note that (1B.2p can be written as 

5 55 = ^2 + -!-)^— + 0(1), 
V u{a) J z - z c 

and if we note that the (5, 5) entry in the Jordan canonical form of the zero order term 
in the Taylor expansion of C(t) is c := C55 = ££55, where t = z — z c , then we can conclude 
that c = (2 + ^y) £ (3, +00) with 

lim c = 3, lim c = +00. 

a— >0+ a— >(a*)- 
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